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ABSTRACT 

We describe the interplay between stellar evolution and dynamical mass loss of evolv- 
ing star clusters, based on the principles of stellar evolution and cluster dynamics and 
on the details of a grid of TV-body simulations of Galactic cluster models. The cluster 
models have different initial masses, different orbits, including elliptical ones, and dif- 
ferent initial density profiles. We use two sets of cluster models: one set of Roche-lobe 
filling models and a new set of cluster models that are initially underfilling their tidal 
radius. 

We identify four distinct mass loss effects: (1) mass loss by stellar evolution, (2) 
loss of stars induced by stellar evolution and (3) relaxation-driven mass loss before 
and (4) after core collapse. At young ages the mass loss is dominated by stellar evo- 
lution, followed by the evolution-induced loss of stars. This evolution-induced mass 
loss is important if a cluster is strongly emersed in the tidal field. Both the evolution- 
induced loss of stars and the relaxation-driven mass loss need time to build up. This 
is described by a delay-function that has a characteristic time scale of a few crossing 
times for Roche-lobe filling clusters and a few half mass relaxation times for ini- 
tially Roche-lobe underfilling clusters. The relaxation-driven mass loss (called "disso- 
lution" in this paper), can be described by a simple power law dependence of the mass 
d(Af/M Q )/d< = — (M/Mq) 1 ~ 1 /to, where to depends on the orbit and environment 
of the cluster. The index 7 is 0.65 for clusters with a King-parameter Wo = 5 for 
the initial density distribution, and 0.80 for more concentrated clusters with Wo = 7. 
For initially Roche-lobe underfilling clusters the dissolution is described by the same 
7 = 0.80, independent of the initial density distribution. The values of the constant to 
are derived for the models and described by simple formulae that depend on the orbit 
of the cluster. The mass loss rate increases by about a factor two at core collapse and 
the mass dependence of the relaxation driven mass loss changes to 7 = 0.70 after core 
collapse. 

We also present a simple recipe for predicting the mass evolution of individual star 
clusters with various metallicities and in different environments, with an accuracy of 
a few percent in most cases. This can be used to predict the mass evolution of cluster 
systems. 
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In this paper we study the mass loss and the mass history 
of star clusters in a tidal field, based on a grid of iV-body 
simulations. The purpose of this study is three-fold: (a) to 
understand and quantitatively describe the different effects 
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that are responsible for mass loss, (b) to study the interplay 
between the different mass loss mechanisms and (c) to de- 
velop a method for predicting the mass history of individual 
clusters of different initial conditions and in different envi- 
ronments. This information is needed if one wants to analyse 
observed star cluster systems in different galaxies. Therefore 
we describe in this paper in detail the different mass loss ef- 
fects; how much each one contributes to the mass loss rate; 
how it depends on cluster parameters and environment and 
how these effects determine the mass history and total life- 
time of a cluster. In particular we will point out the impor- 
tance of the loss of stars that is induced by stellar evolution. 
This mass loss is proportional to the evolutionary mass loss 
and therefore adds to the mass loss rates of clusters at young 
ages. We will show that both the evolution-induced mass loss 
and the relaxation-driven mass loss start slowly with a delay 
time on the order of a few crossing times at the tidal radius. 
We will also show that the mass loss rate after core collapse 
is about a factor two higher than before core collapse, and 
that the dependence of the relaxation-driven mass loss on 
mass is different before and after core collapse. 

The mass of star clusters decreases during their lifetime, 
until they are finally completely dissolved. The stars that are 
lost from clusters add to the population of field stars. The 
mass loss is due to stellar evolution and to several dynami- 
cal effects such as two-body relaxation, tidal stripping and 
shocks. These effects have all been extensively studied indi- 
vidually. However, to understand and describe the combina- 
tion and interplay of the effects one has to rely on dynamical 
simulations. 

The effects of stellar evolution on star clusters can be 
studied by means of stellar evolution tracks for a large 
range of masses and metallicities (e.g. Anders & Fritze-v. 
Alvensleben 2003; Bruzual & Chariot 2003; Fioc & Rocca- 
Volmerange 1997; Leitherer et al. 1999; Maraston 2005). The 
dynamical effects of cluster evolution have been described in 
a large number of theoretical studies starting with Ambart- 
sumian (1938), Spitzer (1940), Chandrasekhar (1943) and 
a series of papers by King, e.g. King (1958). This was fol- 
lowed by the seminal works of Spitzer (1958 and 1987) and 
by many other studies, e.g. Chernoff & Weinberg (1990); 
Gnedin & Ostriker (1997); Aarseth (1999); Fukushige & 
Heggie (2000). The first JV-body simulations of clusters were 
done by von Hoerner (1960). For a review of early iV-body 
simulations of star clusters, see Aarseth & Lecar (1973). 

The recent advancement of computational power, in 
particular the development of the GiMP-E-computers 
(Makino et al. 2003) and the use of Graphics Processing 
Units (GPU) has allowed the improvement and verification 
of these theoretical models by means of direct A r -body sim- 
ulations (Vesperini & Heggie 1997; Portegies Zwart et al. 
1998; Baumgardt & Makino 2003, hereafter BM03; Gieles & 
Baumgardt 2008) . For the purpose of the present study the 
following results are particularly important: 
(i) The realization that mass loss by tidal effects does not 
only scale with the half-mass relaxation time (as was as- 
sumed in earlier studies) , but by a combination of the half- 
mass relaxation time and the crossing time (Fukushige & 
Heggie 2000; Baumgardt 2001; BM03). This implies that 



the lifetime due to evaporation in a tidal field does not scale 
linearly with the cluster mass M, but with A/ 7 with 7 ~ 0.6 
to 0.7. 

(ii) The realization that mass loss by shocks due to the pas- 
sage of spiral arms and giant molecular clouds scales with 
the density M/r 3 of the clusters (Spitzer 1958, Gnedin & 
Ostriker 1997). Adopting the observed mean mass-radius 
relation of r k M 0,1 for clusters in spiral galaxies (Larsen 
2004, Scheepmaker et al. 2007) then also results in a mass 
loss rate that scales approximately as M 1 , with 7 similar to 
the value of evaporation in a tidal field (Gieles et al. 2006, 
2007) . This is in agreement with empirical determinations of 
7 ~ 0.6 from studies of cluster samples in different galaxies 
(Boutloukos & Lamers 2003; Gieles et al. 2005; Gieles 2009). 
(iii) A grid of cluster evolution models with different initial 
masses, different initial concentration factors and in differ- 
ent Galactic orbits by means of iV-body simulations (BM03) 
allows a study of the interplay between stellar evolution and 
dynamical mass loss, that is not easily done by theoretical 
studies. In particular it shows how the mass loss depends on 
mass, age and external conditions, and how the stellar mass 
function evolves during the life of the cluster. 

In this paper we will use a grid of A r -body simulations of 
Roche- lobe filling models (BM03), supplemented with a new 
grid for Roche-lobe underfilling models, of Galactic clusters 
of different initial mass, different initial concentrations and 
in different orbits to describe the process of mass loss from 
clusters and the interplay between the different effects. We 
also derive a method for calculating the mass loss and mass 
history for clusters of different metallicity and in different 
environments. This results in an improvement of the ana- 
lytical description of the mass history of clusters that was 
based on a combination of stellar evolution and dynamical 
effects. (Lamers et al. 2005). 

The paper is arranged as follows. In Sect. 2 we describe 
the mass loss processes of star clusters: stellar evolution and 
dynamical effects. In Sect. 3 we describe the results of N- 
body simulations of BM03 used in this study. Sect. 4 deals 
with the mass loss due to stellar evolution, i.e. both the di- 
rect mass loss and the evolution-induced loss of stars. In 
Sections 5 and 6 we describe the relaxation-driven mass loss 
respectively before and after core collapse. Section 7 deals 
with the mass evolution of clusters in elliptical orbits around 
the galaxy and Sect. 8 deals with initially Roche-lobe under- 
filling clusters. In Sect. 9 we study the relation between the 
total age of a cluster and the initial parameters. In Sect. 10 
we predict the mass loss history of clusters and its main con- 
tributions. Sections 11 and 12 contain a discussion and the 
summary plus conclusions of this study. In Appendix A we 
present a recipe to predict the mass history of star clusters 
in different environments and with different metallicities. In 
Appendix B we tabulate numerical coefficients to calculate 
the mass loss of clusters by stellar evolution. 



2 MASS LOSS PROCESSES 

Clusters lose mass by stellar evolution and by dynamical 
effects, such as two-body relaxation and tidal stripping of 
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stars in a cluster that is emersed in a steady tidal field and 
shocks. The mass loss by stellar evolution is in the form 
of gas ejected by stellar winds and by supernovae, but also 
in the form of compact remnants that may be ejected if 
they get a kick velocity at birth. Mass loss by dynamical 
effects is always in the form of stars. Throughout this paper 
we will refer to these two effects respectively as "mass loss 
by stellar evolution" and "dissolution", either in a steady 
potential field or due to tidal perturbation (shocks). 



2.1 Mass loss by stellar evolution 

The mass fraction that is lost by stellar evolution depends on 
the metallicity and on the adopted stellar initial mass func- 
tion. We have calculated these for clusters with a Kroupa 
IMF, using the evolutionary calculations of Hurley et al. 
(2000) by assuming no dynamical mass loss. The data are 
provided by Pols (2007, Private Communication). 

The various contributions to the evolutionary mass loss 
for (non dissolving) clusters with metallicities of Z=0.0004, 
0.001, 0.004, 0.008 and 0.02 can be expressed with very high 
accuracy (better than ~ 1%) by 3rd order polynomials as 
function of time. We have calculated these fit formulae for 
models with a Kroupa (2001) IMF in the range of 0.1 to 100 
Mq. These models have an initial mean stellar mass of 0.638 
Mq. The fit formulae for clusters are listed in Appendix B 
for the following parameters: 
(i) the remaining mass fraction /i(£) = M(t)/Mi, 
(ii) the mass fractions of black holes [ibh = Mbh/Mi, neu- 
tron stars /iNS and white dwarfs /xwd, 

(iii) the mean mass of all stars < m > and of the black 
holes < m >bh, neutron stars, < m >bns and white dwarf 
< m >wd, 

(iv) the luminosity L(t)/L@ of a cluster with an initial mass 
of 1 M Q . 

The mass fraction that is lost by winds and supernova 
ejecta is 1 — ^(t). If compact remnants are ejected with a 
kick velocity then the remaining mass fraction due to stellar 
evolution is 

H m (t) = fj,(t) - /™ k MBH(t) - /kldcMNs(f) - /k£kMWD(t) (1) 

where f§^ k , /idfk an d fkick are the fractions of these stellar 
remnants that are ejected out of the cluster by their kick 
velocity. If all BHs are kicked out then /j^ = 1 and if 
all WDs are retained then /■$&, = 0. The fraction of the 
luminous mass that is left by stellar evolution is 



f™ m (t) = n{t) - /zbh(£) - Hns{t) - (J,wr>(t). 



(2) 



All fraction /i are expressed relative to the initial cluster 
mass Mi. 

The mass loss rate of a cluster due to stellar evolution 
can now be expressed as 



(dM\ _ „,_ M d Mcv (i) 



V dt A 



Mi um (t) ■ 



dt 



(3) 



which is negative, since (j, BV decreases with time. This expres- 
sion is strictly valid for the early phases of cluster lifetime 
before the preferential loss of low mass stars by dynamical 
effects has changed the shape of the mass function. Since 



stellar evolution dominates the mass loss only in the early 
phase of the clusters lifetime, equation [3] is a good approx- 
imation. (For a description of evolutionary mass loss in a 
cluster with preferential loss of low mass stars see Kruijssen 
& Lamers 2008, Kruijssen 2009 and Trenti et al. 2010.) 

In the description of the mass loss of the cluster models 
studied by iV-body simulations in this paper the effect of the 
changing mass function due to evolution and the preferential 
loss of low mass stars is properly taken into account, as it is 
in the output of the simulations. 



2.2 Mass loss by dynamical effects or 
"dissolution" 

The time-dependent mass loss by dissolution can be 
described by (dM/dt)dis = —M/tdis, where tdis = 
(dln(M)/dt) is the dissolution time scale that depends on 
the actual cluster mass and on the environment of the clus- 
ter. Let us assume that we can describe tdis as a power-law 
function of mass, as tdis = to(M / Mq) 1 , with the constant to 
being the dissolution parameter (which is the hypothetical 
dissolution time scale of a cluster of 1 Mq). The changes in 
the cluster mass due to dissolution is then described bj£j 



fdM\ 
V dt Jdis 



M(t) 

'tdis(t) 



Mjt) 1 - 



(4) 



We stress that the dissolution time scale tdis is not the same 
as the total life time of the cluster, ttot, although these are 
related. Integration of Eq. [4] shows that ttot = toM 1 /"/ in the 
absence of stellar evolution, where Mi is the initial mass. In 
reality stellar evolution also removes part of the cluster mass. 
This implies that this simple estimate of ttot overestimates 
the real cluster lifetime. 

Theoretical considerations suggest that 7 ~ 0.65 to 
0.85. This follows from the following dynamical arguments. 



2.2.1 Dissolution in a steady tidal field 

Spitzer (1987) has argued that a fraction of the stars £ es- 
capes each t r h, such that 

(dN\ _ jN 

V dt / dis t r h 



(5) 



where t rh oc (7V 1/2 /ln A)r^ i/2 m~ 1/2 is the half-mass relax- 
ation time, rh is the half-mass radius, m is the mean stellar 
mass, N — M/m is the number of stars and In A is the 
Coulomb logarithm. The value of £ is larger for Roche-lobe 
filling clusters than for clusters in isolation. 

In analytical studies of cluster dissolution a single value 
for £ is usually assumed, implying that the cluster life- 
time is a constant times t r h (see e.g. Spitzer 1987). How- 
ever, a recent theoretical study by Lee (2002) and Gieles 
& Baumgardt (2008) has shown that £ is a strong function 
of the Roche- lobe filling factor rh/rj, (where rj is the Ja- 
cobi radius, i.e. the Roche-lobe radius for clusters). These 

1 Throughout the rest of this paper all masses M; , M and m are 
in units of Mq and all ages are in Myrs. 
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authors found for clusters with r^/rj ^ 0.05 that £ scales 



roughly as (rh/rj) ' . Since rj oc M ' w~ , where uj is 
the angular frequency of the cluster orbit, f oc t C ross x w. 
Here t cr oss is the mean crossing time of stars in a cluster: 
icross oc rh /vGM .This dependence of £ on r^ 3 ' 2 cancels 
the rh 3 ' 2 dependence of t r h such that the radius becomes an 
unimportant parameter in AN /At. This can be understood 
intuitively as follows: for a smaller (larger) radius, relaxation 
becomes more (less) important, while the escape criterion 
due to the tidal field becomes less (more) important. 

On top of this, Baumgardt (2001) and BM03 showed 
that the dissolution time scale does not scale linearly with 
t r h, but rather with a combination of t r h and the crossing 
time, tcross, because even unbound stars need time to leave 
a cluster (see Fukushige & Heggie 2000 for details). 

The relevant time scale is 



^dyn OC £ r h £ cross 

N 



WaJ 



(6) 
(7) 



This implies that dissolution time scale in terms of number 
of stars is 



idis 



N 



AN /At 

N 



^dyn 



VlnAj 



(8) 
(9) 



where idis^ is the dissolution time scale if dissolution is 
expressed in terms of N instead of M. In the range of 
10 4 < N < 10 6 we can approximate A ~ Q.02N (Giersz 



& Heggie 1994) and JV/lnA oc AT ' 80 and so tdis 



N p 



with p = 0.80x. 

In this paper we describe the dissolution time as a func- 
tion of M instead of N. Comparing the two expressions 
tdis oc M 7 and tdis^ oc N p we see that 7 7^ p if the mean 
stellar mass m changes during the clusters lifetime. If m de- 
creases as function of time, i.e. increases as function of M(t), 
then 7 < p, whereas 7 > p if m increases with time. We will 
see below that after an initial phase, dominated by stellar 
evolution, m increases with time due to the preferential loss 
of low mass star by tidal stripping. So the values of 7 are 
expected to be slightly higher than the values of p. 

2.2.2 Dissolution due to shocks 

Clusters can also be destroyed by shocks (e.g. Ostriker et al. 
1972, Spitzer 1987, Chernoff & Weinberg 1990, Gnedin & 
Ostriker 1997) due to encounters with spiral arms or giant 
molecular clouds in the disk of a galaxy, by disk shocking 
for clusters in orbits inclined with respect to the Galactic 

2 Using the total lifetime as an indicator of the dynamical time 
BM03 found that x = 0.75 for Roche-lobe filling models with an 
initial concentration factor of the density King-profile Wo = 5 
and x = 0.82 for the more centrally concentrated Wo = 7 models. 
This would imply p ~ 0.60 and 0.66 for Roche-lobe filling models 
of Wo = 5 and 7 respectively. 



plane and by bulge shocking for clusters in highly elliptical 
orbits. 

The mass loss rate due to shocks depends on the cluster 
half-mass radius, rh, as (AM/ At) oc M/ ' p\ oc rh 3 , with ph the 
density within rh. This implies that the mass loss time scale 
depends on the cluster properties as tdis = —M/(AM /At) oc 
M/rh 3 , which is proportional to the cluster density. The 
observed mean mass-radius relation of clusters is not well 
defined, but Larsen (2004) and Scheepmaker et al. (2007) 
find a mean relation of rh oc M x with A ~ 0.13. So the 
time scale for mass loss due to shocks is tdis oc M ' 61 for a 
constant mean stellar mass (Gieles et al. 2006, 2007). This 
dependence is almost the same as that for tidal dissolution. 



2. 2. 3 The expected values of to and 7 

Based on the arguments of the previous subsections we ex- 
pect that the combined mass loss by tidal dissolution and 
shocks can be described by a function of the form 



V dt Jdi, 



N 

'tdis* 



lo^ 



V dt /dis 



Al 



f 1 — r 



to 



-(10) 



with p ~ 0.65 and 7 <; p. We will use the second expression 
with the value of 7 derived from the A r -body simulations of 
BM03. 

In an environment where cluster dissolution is only 
due to stellar evolution, internal dynamical effects and tidal 
stripping, to depends on the potential field in which the clus- 
ter moves. If clusters move in elliptical orbits in a logarith- 
mic potential field, i.e. a constant galactic rotation veloc- 
ity u Gai, the dissolution time is reduced by a factor 1 — e, 
compared to a circular orbit at Ra, where the eccentricity 
e = (Ra — Rp)/{Ra + Rp) and Ra and Rp are the apogalac- 
tic and perigalactic distances respectively. This implies that 
we expect the value of to to vary as 



to 



(i—i) ( flGal ^ ( 
V WP J ^8.5kpcy V 



"Gal 



220km s" 



(11) 



where t rcf is a constant, whose value will be derived from 
the A^-body simulations of BM03. 

The factors 7?Gai/8.5kpc and iiGai/220km s" 1 provide 
a scaling for calculating the dissolution of clusters at differ- 
ent galactocentric distances in other galaxies with constant 
rotation velocity. The factor 1/m 7 is a result of the conver- 
sion of Af to M if the mean stellar mass m is about constant. 
We will see below that this yields a very good description of 
the mass loss rate and M(t) for all models of BM03, if m is 
chosen appropriately!^ 



3 If other processes, such as encounters with GMCs or spiral den- 
sity waves, are important, then the value of to will be smaller than 
predicted by Eq. [TTJ Wielen (1985) and Lamers & Gieles (2006) 
found that clusters in the solar neighbourhood are mainly de- 
stroyed by encounters with GMCs which reduces to by a factor 4 
compared to Eq. 1111 Kruijssen & Mieske (2009) have derived the 
values of to for a number of galactic globular clusters in elliptical 
orbits, assuming 7 = 0.70. 
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Table 1. The N-body models of BM03 used in this study 



Input parameters 












Timescales 








Output parameters 












# 


Mass 


nr 


Wo 


flGal 


Orbit 


n 


n. 


h% 


tec 


trh 


*cr (r h ) 


tcr(n) 


7 


to 


t M? 


tdel 


fsr 


*rj C 


jee 




M 


stars 




kpc 




pc 


pc 


Gyr 


Gyr 


Gyr 


Myr 


Myr 




Myr 


Gyr 


Myr 




Myr 




i 


71952 


128k 


5 


15 


circ 


89.6 


16.75 


45.3 


36.42 


7.20 


15.24 


133 


0.65 


40.0 


57.4 


400 


4.1 


16.1 


1.60 


2 


35915 


64k 


5 


15 


circ 


71.0 


13.28 


26.9 


22.51 


3.88 


15.24 


133 


0.65 


42.0 


38.4 


400 


1.0 


17.2 


1.01 


3 


18205 


32k 


5 


15 


circ 


56.7 


10.59 


19.8 


14.24 


2.13 


15.24 


133 


0.65 


41.2 


21.2 


400 


1.0 


20.0 


1.40 


4 


8808 


16k 


5 


15 


circ 


44.5 


8.32 


13.4 


8.72 


1.13 


15.24 


133 


0.65 


37.9 


13.7 


400 


0.7 


20.4 


1.29 


5 


4489 


8k 


5 


15 


circ 


35.5 


6.64 


9.0 


5.92 


0.63 


15.24 


133 


0.65 


36.0 


8.5 


100 


0.7 


19.0 


1.37 


6 


71236 


128k 


5 


8.5 


circ 


61.1 


11.43 


26.5 


21.34 


4.05 


8.63 


76 


0.65 


21.5 


32.1 


228 


0.8 


8.5 


1.70 


7 


36334 


64k 


5 


8.5 


circ 


48.8 


9.13 


17.2 


13.91 


2.22 


8.63 


76 


0.65 


22.0 


20.3 


228 


0.8 


10.4 


1.40 


8 


18408 


32k 


5 


8.5 


circ 


39.0 


7.28 


11.1 


8.41 


1.22 


8.63 


76 


0.65 


21.7 


12.8 


228 


0.8 


9.9 


1.50 


9 


9003 


16k 


5 


8.5 


circ 


30.7 


5.74 


7.5 


5.06 


0.65 


8.63 


76 


0.65 


20.5 


7.4 


228 


0.7 


10.7 


1.34 


10 


4497 


8k 


5 


8.5 


circ 


24.3 


4.55 


4.9 


3.30 


0.36 


8.63 


76 


0.65 


20.0 


4.7 


228 


0.8 


10.2 


1.42 


11 


71218 


128k 


5 


2.8 


circ 


29.4 


5.50 


9.3 


7.66 


1.35 


2.87 


25 


0.65 


7.5 


10.7 


75 


0.7 


3.0 


1.02 


12 


35863 


64k 


5 


2.8 


circ 


23.4 


4.37 


5.9 


4.63 


0.73 


2.87 


25 


0.65 


6.7 


5.8 


75 


0.7 


3.3 


1.36 


13 


18274 


32k 


5 


2.8 


circ 


18.7 


3.49 


3.6 


2.85 


0.40 


2.87 


25 


0.65 


0.0 


3.5 


75 


0.6 


3.3 


1.26 


14 


9024 


16k 


5 


2.8 


circ 


14.8 


2.76 


2.3 


1.58 


0.22 


2.87 


25 


0.65 


5.3 


2.0 


75 


0.4 


3.2 


1.16 


15 


4442 


8k 


5 


2.8 


circ 


11.7 


2.18 


1.3 


0.85 


0.12 


2.87 


25 


0.65 


4.4 


1.0 


75 


0.4 


2.8 


1.13 


16 


71699 


128k 


7 


8.5 


circ 


61.3 


7.11 


28.5 


12.62 


1.99 


4.22 


76 


0.80 


6.4 


40.1 


228 


0.8 


11.5 


1.52 


17 


35611 


64k 


7 


8.5 


circ 


48.5 


5.63 


17.2 


7.87 


1.07 


1.22 


76 


0.80 


6.5 


28.2 


22X 


0.8 


10.5 


1.57 


18 


18013 


32k 


7 


8.5 


circ 


38.7 


4.48 


11.2 


4.87 


0.58 


1.22 


76 


0.80 


6.5 


15.7 


228 


0.8 


10.5 


1.43 


19 


8928 


16k 


7 


8.5 


circ 


30.6 


3.55 


6.9 


2.89 


0.32 


1.22 


76 


0.80 


0.0 


8.5 


228 


0.8 


11.0 


1.21 


20 


4402 


8k 


7 


8.5 


circ 


24.2 


2.80 


4.4 


1.67 


0.17 


1.22 


76 


0.80 


5.5 


4.3 


228 


0.8 


10.1 


1.14 
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Left section: model parameters; Middle section: cluster time scales; Right section: fit parameters. 
First three blocks: clusters with Wo = 5 in circular orbits at i?Gal= 15, 8.5 and 2.8 kpc. Fourth block: clusters with Wo = 7 in circular 

orbits at i?Gal=8.5 kpc. Fifth block: clusters with Wo = 5 in elliptical orbits with apogalactic distance of Ra = 8.5 kpc and a 

perigalactic distance of Ra{^ — e)/(l + e )- The values of rh and rt apply to the perigalacticon. Clusters with Wo = 5 or 7 have 7 = 0.65 

and 0.80 respectively, before core collapse. The number of stars is given in units of lfc = 1024. 



3 THE MODELS OF ROCHE-LOBE FILLING 
CLUSTERS 

In this paper we study the results of iV-body simulations, in 
order to understand the way clusters evolve due to dynam- 
ical and evolutionary effects. We use the models of Roche- 
lobe filling clusters from BM03 for comparison with our pre- 
dictions, supplemented with a few models of initially Roche- 
lobe underfilling clusters. Out of the initial 33 models we 
have selected 25 representative BM03 cluster models. These 
are chosen because they allow the study of the effects of 
the different parameters of the models, i.e. initial mass and 
initial concentration, and of the cluster orbits, i.e. Galacto- 
centric distance and eccentricity. (We found that the infor- 
mation derived from these models also applies to the other 
8 models.) The models span a range of dissolution times 
between 1.5 and 50 Gyr. The selected models are listed in 
Table [TJ 



3.1 Parameters of the cluster models 

The models are divided into five blocks, separated by hori- 
zontal lines in Table [l] The first three blocks contain mod- 
els of Roche-lobe filling clusters of different masses with an 
initial concentration factor Wo = 5, in circular orbits at 
galactocentric distances of Rami = 15, 8.5 and 2.83 kpc. The 
fourth block contains Roche-lobe filling cluster models in 
circular orbits and with different initial masses, but with a 



more concentrated initial density distribution, with a King 
profile of Wo = 7. These will be used to study the effect of 
the initial concentration on the cluster evolution, by com- 
paring them with the results of the Wo = 5 models. The 
fifth block contains Roche-lobe filling cluster models with 
Wo — 5 in elliptical orbits with various eccentricities. All 
models have a Kroupa stellar initial mass function (IMF), 
with a mass range of 0.15 to 15 M@, an initial mean stel- 
lar mass of 0.547 M e and a metallicity of Z = 0.001. The 
clusters have no initial binaries, but binaries do form during 
the dynamical evolution, mainly in the high density central 
region during core collapse. Neutron stars and white dwarfs 
are retained in the cluster when they are formed (no kick 
velocities) but may be lost later by dynamical effects. Black 
holes are not considered in the BM03 models, because of the 
adopted upper mass limit of 15 Mq. 

The clusters in elliptical orbits (nrs 21 to 25) have about 
the same initial mass and the same apogalactic distance of 
Ra — 8.5 kpc, but different elliptical orbits with eccentric- 
ities 0.2 ^ e ^ 0.8. This implies perigalactic distances of 
Rp = Ra(1 - e)/(l + e) between 0.667 Ra and 0.111 Ra. 
For these models the mass loss rates are strongly variable 
with time: at perigalacticon the rates are much higher than 
at apogalacticon. The initial values of the tidal radius, half- 
mass radius etc. in Table [l] refer to the values at perigalac- 
ticon. 

The data of each model in Table [l] are given in three 
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groups, separated by a vertical line. The left group gives 
the initial model data: model nr, initial mass, initial nr of 
stars, Wo, apogalactic distance, type of orbit with eccen- 
tricity, tidal radius r t and half mass radius r a . The middle 
group gives the various time scales of the models: the time 
t\% when M(t) = O.OlAfj, the core-collapse time t cc , the 
initial half mass relaxation time £ rn , the half mass crossing 
time i cr ('"h) and the initial crossing time at the tidal radius 
icr(r-t). 

The right hand group gives the data that describe the 
mass loss rates of the models: to, toM^ (which is a proxy of 
the expected total life time). The values of delay-time tdciay 
and /iJ'd' 1 together describe how the clusters react dynami- 
cally to mass loss by stellar evolution (see Sect. [6]). The last 
two columns give the values of to c , which describes the mass 
loss rate after core collapse, and j cc which describes the in- 
crease in mass loss due to core collapse. The determination 
of these parameters is described below. 

3.2 The mass loss rates of the cluster models 

BM03 define a star to be lost from a cluster if it is outside 
the Jacobi radius rj of the cluster. Stars with a velocity 
v > v csc but r < r.j are still considered cluster members. 
On the other hand, the mass lost by stellar evolution (i.e. 
by winds and supernovae) is assumed to leave the cluster 
immediately. 

We have derived the mass loss rates of the N-body mod- 
els of BM03. From the output of these model calculations 
we can separate the mass loss that is due to stellar evolu- 
tion from the contribution by dynamical effects. There is 
an important difference between these two mass loss rates. 
Mass loss by stellar evolution is instantaneous and indepen- 
dent of the structure and orbit of the cluster. On the other 
hand, mass loss by dynamical effects always proceeds on a 
slow time scale and needs time to build up. We will see this 
in the results. 



3.3 Three phases of mass loss 

A study of the mass loss rates of the BM03 models shows 
that three mass loss phases can be recognized. This is de- 
picted in Fig. [T]for two models (nrs 6, 12), which shows the 
variation of {AM /At) as function of M and t. In the first 
phase (A) mass loss is dominated by stellar evolution and 
so the mass loss rate drops steeply with time. In the sec- 
ond phase (B) mass loss is dominated by dynamical effects 
and the mass loss rate behaves approximately as a power 
law of {AM /At) oc M 1-7 (see Sects. [23T1 and [2T2~2|) . The 
third phase (C) is after core collapse. Mass loss is also dom- 
inated by dynamical effects, but the mass loss rate is higher 
than before core collapse. This has been noticed before in 
the mass evolution of cluster models, e.g. Baumgardt (2001). 
The separation between the three regions is not as strict as 
suggested in Fig. \T\ because the different effects overlap near 
the boundaries (see below). 

In the next three sections we discuss the three mass loss 
phases and how they depend on the cluster parameters and 
the environment. 



4 DIRECT AND INDUCED MASS LOSS BY 
STELLAR EVOLUTION 

Mass loss during the early history of clusters (phase A in 
Fig. [l]) is dominated stellar evolution. 

At very early stages only stellar evolution contributes to 
the mass loss. But very soon thereafter, the mass loss rate 
{AM /At) of all models is higher than (dA//di) cv . This is 
seen best in the right hand panels of Fig. [TJ where the mass 
loss rates are plotted versus time: the difference between 
the total mass loss rate (top lines) and the mass loss by 
evolution (dotted line) is much larger than the mass loss by 
dissolution (almost straight line), which will be discussed 
below. This shows that during the early phases, when the 
mass loss is dominated by stellar evolution, Roche-lobe filling 
clusters lose an extra amount of mass (in the form of stars) 
by dynamical effects, induced by the mass loss due to stellar 
evolution. This evolution-induced mass loss is due to the fact 
that the cluster radius expands and the tidal radius shrinks 
due to evolutionary mass loss. 

For adiabatic models, we expect that the evolution- 
induced mass loss rate, {AM/At)^ d will be about equal to 
the mass loss rate by stellar evolution, (dM/dt) cv . This can 
be understood as follows. 

If a cluster loses a fraction S « 1 of its mass Mq 
on a time scale longer than its crossing time, the clus- 
ter will adiabatically expand such that its radius, relative 
to its initial radius, is r/ro = Mo/M ~ 1 + 5. At the 
same time, the mass loss causes the Jacobi radius to shrink: 
rj/rjo = (M/M ) 1/3 w 1 - 8/3. The mass in the shell be- 
tween 1 — 8/3 < r/rjo < 1 + 8 is consequently unbound 
since it is outside the new Jacobi radius. For a logarithmic 
potential, the density at rjo, p{rjo), is 6 times lower than 
the mean density within rjo: pjo — 3Mo/{Aitr 3 JQ ). The mass 
AM that is in the unbound shell is thus AM = A-Kr 2 J0 pjoAr, 
with Ar = ASrjo/3, such that AM = (2/3)5M . In fact, 
AM will be slightly higher because we have adopted the 
lower limit for the density in the outer layers of the clus- 
ter. So we may expect that evolutionary mass loss rate in- 
duces about the same rate of evolution-induced mass loss, 
(dM/dt)& = find X (dM/dt)ev with / ind ~ 1. 

A study of all models shows two deviations from this 
simple expectation. 

(i) The value of /i n d is only about unity for models for which 
(dM/dt) ev >> (dAf/dt)di s - This is the case for models with 
a very long lifetime of itot > 25 Gyr. If the mass loss by disso- 
lution in the early lifetime of the cluster cannot be ignored, 
the evolution-induced mass loss rate is smaller. Therefore 
/ind will be smaller than unity. 

(ii) The evolution-induced mass loss does not start at t = 
but it needs time to build up. We can expect that the time 
scale for this build-up will be of the order of the crossing 
time at the tidal radius, because this is the time scale on 
which stars can leave the cluster by passing the tidal radius 
due to the reduction of the depth of the potential well. 

A study of all the cluster models shows that we can 
describe the evolution-induced mass loss rate by 

(^): d - A-W x (^) cv . (12) 

with 
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Figure 1. The three phases of mass loss in two Af-body models of BM03: A = dominated by stellar evolution, B = dominated by 
dissolution. C = dominated by dissolution after core collapse. The line separating phases A and B is when the evolutionary mass loss 
has dropped to 10% of the total mass loss. The line separating phases B and C indicates the core collapse time. The figure shows dM/dt 
versus M (left) and t (right) in logarithmic units. The full upper line is the total mass loss rate; the dashed line is the mass loss by 
stellar evolution; the full smooth lines are the mass loss by dynamical effects, assumed to scale as a power law of M, before and after 
core collapse. The models are defined by a vector containing: model nr, nr of stars, total lifetime (in Gyrs), concentration parameter Wo, 
Rg (in kpc), and orbit. 



/ind(t) = /i„d X X /delay (t) 



(13) 



We assume that the growth of the evolution-induced mass 
loss approaches its maximum value in an exponential func- 
tion of the delay time scale 



/delay (£) = 1 — exp{ — (t/t delay ) } 



(14) 



The delay time scale is expected to depend on the crossing 
time at the tidal radius, so 



^delay <Melay " ^cr\' t) 



(15) 



Eq.[l3]describes the increase of fi n d(t) from to /^d" with a 
time scale of tdelay Figure [2] shows the function /i n< j for two 
models. The figure shows that the exponential expression 
describes the delay function quite well. 

The values of Wdciay turn out to be about 3.0 for all 
Roche-lobe filling models. Therefore we have adopted this 
value for all models, (see Table [TJ. 

The values of fg£* listed in Table[[]show that /££* = 1 
only for the clusters with very long lifetimes larger than 
about 25 Gyr. For clusters with shorter lifetime /j^ d x does 
not reach this value, because the contribution by dissolu- 
tion helps to restore the equilibrium that was destroyed by 
the fast mass loss due to stellar evolution. Therefore we ex- 
pect the value of /i^d X to depend on the ratio between the 
evolutionary mass loss, (dM/di) ov = M(d/x/dt)ev and the 
mass loss by dissolution, —M/tdis with idis = toM 1 . The 
smaller the ratio (dM/di) dis /(dM/dt)*, the larger /™ 
with a maximum of /™d* — 1- Fig. [3] shows the values of 



/hid X as a function of toM^ . (The determination of the val- 
ues of to and 7 is described in Sect. 5). The results can be 
represented by the relation 



./ ii; 



-0.86 + 0.40 xlog(toM^) for ]og(t M?) ^ 2.15(16) 



and /™d x = if log(toM^) < 2.15, with the ages in Myrs. 



5 DISSOLUTION BEFORE CORE COLLAPSE 

We first consider the dynamical mass loss in phase B, i.e. be- 
fore core collapse, which is expected to vary as (dM/dt)dw ~ 
M 1 ' (Sect. [2~2)l . This mass loss is due to two-body relax- 
ation (dissolution) plus the induced loss of stars due to the 
expansion of the cluster and the shrinking of the tidal radius. 
In order to derive the value of 7 we study the dependence 
of (dM/di)di s on M (t) for all models of clusters in circular 
orbits. For the determination of the mass loss by dissolu- 
tion we have corrected the total mass loss rate for both the 
evolutionary and the evolution-induced mass loss rate. 



5.1 The empirical value of 7 before core collapse 

The upper left panel of Fig. [4] shows the mass loss rates as 
a function of mass of models 1 to 15, i.e. for clusters with 
Wo — 5. For each model the mass loss rates are measured in 
the mass interval where the evolutionary mass loss is small, 
(dM/dt) cv < 0.1(dA//di), up to core collapse. So each model 
occupies a limited region in the M(t)-range. The mass loss 
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Figure 2. The function /delay for two N-body models with very 
different delay time scales. The smooth line is the fit by formula 
1141 The parameters are given in the figure. 
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Figure 3. The values of / m J x versus logtgM? for all Roche-lobe 
filling cluster models. For models in circular orbits the values of 
J™ decrease with decreasing t$M? , i.e. with increasing mass 
loss by dissolution, as expected. The dashed line shows the ap- 
proximate relation which is given by Eq. 1161 Clusters in highly 
eccentric orbits have smaller values of /^^"(see Sect. 0. 



rates are corrected for evolutionary and evolution-induced 
mass loss. 

The mass loss rates are normalized to RgbX = 8.5 kpc 
and to the mean stellar mass, to make the curves overlap. 
This is possible because we expect from Eqs. llOl and llll that 



(dM\ /_Rcai_ 
V dt / I 8.5kpc 



(£)-(£)«"■ 



(17) 



where m is the mean stellar mass and Wi is the initial 
mean stellar mass. We see that the curves nicely overlap. 
We have fitted a straight line through the data with a slope 
1 - 7 = 0.35 ± 0.02. This implies that 7 = 0.65 ± 0.02. (The 
appearance of 7 in the term (jn/mi) 1 ~ 1 implies that we 
had to derive the value of 7 in two iterations.) We will adopt 
7 = 0.65 for the Wo = 5 models in the rest of the paper. 
This agrees with observations of star clusters in M51 (Gieles 
2009). and with the N-body simulations of clusters without 
stellar evolution by Gieles & Baumgardt (2008). 

The middle left panel shows the relation for the nor- 
malized mass loss rates for the models with Wo = 7. Each 
model occupies only a small part in this diagram because the 
pre-core collapse time of these models is short: only about 
half as long as that of the Wo = 5 models. Nevertheless, we 
see that there is a clear power law dependence on M with a 
slope of 1 - 7 = 0.20 ± 0.04 which implies 7 = 0.80 ± 0.04. 
We adopt 7 = 0.80 for the pre core collapse phase of the 
Wo = 7 models. This value is higher than for the Wo — 5 
models, as expected (see Sect. I2.2~T]| . 

5.2 The dissolution parameters to and i^ f 

We have fitted the mass loss rates of the BM03 models in 
phase B, i.e. the pre-core collapse phase, with a function 
(dM/dt)dis = -M 1 ~ 7 /t with 7 = 0.65 and 0.80 for the 
models with Wo — 5 and 7 respectively. The values of the 
mean stellar mass in that phase is also derived from the 
details of the BM03 models. The values of to are listed in 
Table [T] We see that they are approximately constant for 
the Wo = 5 cluster models at the same value of i?Gai and 
that for clusters in different orbits to oc i?Gai, as expected 

(Eq.rnj. 

Fig. [5] (top panel) shows the relation between to and 
toM 7 , which is a proxy for the total lifetime of the cluster, 
for all models. For each value of i?Gai the values of to de- 
crease slightly with decreasing to A/ 7 . This is because to is 
expected to depend on m -7 (Eq. II 1ft and the mean mass in 
the pre-core collapse phase depends on toM 7 . The vertical 
offset between the relations for different values of i?Gai is 
because to is also expected to be proportional to i?Gai (Eq. 
Hip . The difference in the values of to between models of 
Wo = 5 and 7 at the same RguI is due to the difference in 
the value of 7. In fact, for clusters with Mi = 1O 4 M0 the 
dissolution time scale tdi s = tolO 47 of a Wo = 5 cluster is 
about the same as that for a Wo — 7 cluster. 

The values of to are measured after the stellar evolution- 
dominated phase but before core collapse. This is when 
M(t)/M[ = /i is approximately between 0.4 and 0.15, with 
a mean value of < /1 >~ 0.25. The middle panel of Fig. [5] 
shows the values of m in the pre-collapse phase at fi = 0.25 
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Figure 4. The dependence of the mass loss rate by dissolution on M(t). The figure shows the values of log — (dM/di)Jj? rm , defined by 
Eg. 1171 versus log(M), before (left) and after (right) core collapse for the models with Wo = 5 (top), Wo = 7 (middle) and Roche-lobe 
underfilling models (bottom). Notice the power law dependence of (dAf/dt)^? rm on M. The dashed lines show the best-fit linear relation. 
Before core collapse these relations have a slope of 1 — 7 = 0.35, 0.20 and 0.20 respectively for Wo = 5, Wo = 7 and Roche-lobe underfilling 
models. After core collapse the relation is expressed by a double power law with 1 — 7 CC = 0.30 if M > 10 3 Mq and 1 — 7 CC = 0.60 if 
M < 10 Mq. The dotted line in the lower right hand panel shows the shape of the expected relation with the Coulomb logarithm taken 
into account if the mean stellar mass is m = 0.5 Mq (see Sect. [6721 . 



derived from the BM03 models. These values of m are ap- 
proximately 



logm 



+0.184 - 0.121 x log t M°- b5 for W = 5 



logm = +0.090 - 0.094 x log to Mi for W = 7 (18) 

These two relations are shown in the figure. The increase of 
m towards shorter lifetimes is because at a given fraction of 
its lifetime the maximum mass of a star in a cluster with a 
short lifetime is higher than in case of a longer lifetime due 
to stellar evolution (see also Sect. 15.31 . 

The lower panel of Fig. [5] shows the value of 
tom 7 (7?Gai/8.5)~ 1 (wGai/220)(l — e) _1 , which is expected to 
be constant, t ref (see Eq. Hip . We see that the values are 
indeed about constant with log(£^ f /Myr) = 1.125 ± 0.016 
for Wo = 5 models and 0.550 + 0.015 for W = 7 models. For 
clusters with total ages less than about 3 Gyrs, the values 
are slightly smaller. This is due to the fact that these clus- 
ters still contain massive stars over most of their lifetime 
and massive stars are effective in kicking out lower mass 
stars (BM03). 



orbits. The lower panel shows the mean mass as function of 
time in terms of t/ti%. The mean stellar mass initially de- 
creases due to the loss of high mass stars by stellar evolution, 
but then increases again with age due to the preferential loss 
of low mass stars by dynamical effects. The minimum value 
is reached for all models at t ~ 0.2^%. This is the time when 
the cluster is fully mass segregated after which the low mass 
stars in the outer part of the cluster are lost preferentially 
(BM03). For models with a short total dissolution time the 
minimum mass is higher than for models with a long life- 
time. Notice that after mass segregation m increases as a 
power-law of \i with m oc ^T b with S ~ 0.30 for all BM03 
models, including the ones not shown here. 

The middle panel of Fig. [5] shows that the mean stellar 
mass of the Wo = 7 models is slightly larger than for the less 
concentrated Wo — 5 models. This is due to the fact that 
core collapse occurs earlier in the more concentrated Wo = 7 
models and so the value of to in the pre-core collapse phase 
refers to an earlier time than in the Wo — 5 models. Fig. 
[6] shows that m increases with age, so it is smaller in the 
pre-core collapse phase of the Wo — 7 models. 



5.3 The evolution of the mean stellar mass 

The time variation of m for the different models is shown in 
Fig. [5] The top panel shows the time evolution of m (plotted 
in terms of n) for all models with Wo = 5 and 7 in circular 



5.4 The start of the dissolution process 

We have derived the mass loss by dissolution for the BM03 
models from the output files of these models. The data show 
that the dissolution does not start at t — but needs time 
to develop, just like the evolution-induced mass loss needs 
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Figure 5. The parameters of the models 1 to 20 of clusters in circular orbits before core collapse (left) and after core collapse (right). 
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time to get going. This can be seen in Fig. [T] which shows 
that at very early times the total mass loss rate is equal 
to the evolutionary mass loss, without any contribution by 
dynamical effects. The time needed to develop the dynamical 
dissolution depends on the crossing time at the tidal radius 
and is expected to behave in the same way as the growth 
of the evolution-induced mass loss rate. Therefore we can 
describe the dissolution before core collapse as 



V&t) = ~ /dolay x ~^T 



(19) 



with /delay given by Eq. [14] 

The delay-time of the dissolution is a result of the initial 
conditions of the cluster models which start without stars 
with escape velocities. This might not be realistic. If clusters 
form in a collapsing cloud and go through a phase of violent 
relaxation, the tail of the Maxwellian velocity distribution 
is initially filled and dissolution will start immediately. 



6 DISSOLUTION AFTER CORE COLLAPSE 

We now consider the dissolution in phase C, i.e. after core 
collapse. Fig. [l] shows that the mass loss rate increases by 
a small factor at about t cc and that the slope 1 — 7 of the 
log— (dM/dt)dis vs logAf relation is different from that in 
phase B. 



6.1 The time of core collapse, i cc 

The core collapse times, listed in Table [l] are derived em- 
pirically from the iV-body calculations. Theory predicts that 
core collapse occurs after about a fixed number of central re- 
laxation times, £ r c, (e.g. Spitzer 1987, p. 95). However as the 
cluster evolves, it loses mass and expands, so t TC changes con- 
tinuously. Therefore we search for an empirical expression of 
tec in terms of the initial values of £ r h and 7?Gai (see Fig. [7]). 
A linear regression analysis shows that we can approximate 
tec of all Roche-lobe filling cluster models, including those 
with Wo = 7 and elliptical orbits, to an accuracy better than 
about 0.05 dex by the following relations 

log (t cc ) = 1.228 + 0.872 x log(i rh ) (20) 



6.2 The value of 7 after core collapse 

Core collapse changes the density distribution of the stars 
in the clusters, so that it becomes independent of the initial 
distribution. Therefore we expect the values of 7 in phase C 
to be the same for all the models. The right hand panels of 
Fig.Ushow the normalized dissolution rate, given by Eq. 1171 
after core collapse for cluster models with Wo = 5 (top) and 
Wo = 7 (middle). Each model occupies a certain mass range, 
starting at core collapse down to about M(t) = 10 2 M©. 

Notice that, apart from a small vertical offset, the sets 
of models show very similar lines with a slight curvature, 
in the sense that the slope gets steeper towards lower mass. 
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Figure 7. The relation between the initial half mass relaxation 
time, t r h and the core collapse time t cc for Roche-lobe filling clus- 
ters, including those in elliptical orbits. 



(see Sect. 12. 2. lj) . The variation of log (dM/di) dis with log M, 
predicted with this aproximation is shown in the lower part 
of Fig. |4]by a dotted line with an arbitrary vertical offset. 
The shape of this predicted line is very similar to the one 
derived empirically. 



Figure 6. Top: The mean stellar mass of Wo = 5 and Wo = 7 
clusters as a function of the remaining mass fraction for models 1 
to 20. The shorter the lifetime of a cluster, the higher the line in 
this figure. The vertical dotted lines indicate approximately the 
range where to was determined in the pre-core collapse phase, 
with the central value indicated by a full line. The diamonds 
indicate the minimum values. The mean mass at core collapse 
is indicated by an asterisk or a cross for models with Wo = 5 
and Wo = 7 respectively. Bottom: The mean stellar mass as a 
function of the fraction of the total cluster lifetime. For all models 
the minimum is reached at t ~ 0.15^%. 



The empirical relations can be fitted very well with a broken 
power-law relation 



(AM\ 
\~~aTJ 



postcc 



M{tf 



-7cc 



(21) 



with 7 CC = 0.70 for M(t) > 10 3 and j cc2 = 0.40, where 
the subscript 2 refers to the mass loss at M(t) < 1O 3 M0. 
These values are the same for Wo — 5 and 7 clusters, be- 
cause core-collapse results in a redistribution of the density 
profile which becomes nearly independent of that in the pre- 
collapse phase. Continuity of the mass loss rate at M{i) = 
10 3 M o requires that t c c2 = t c c x I0 3 <7cc-7cc2) = 1Q o.90 ^cc 

The steepening of the slopes in Fig.[4]is the result of the 
changes in the Coulomb logarithm towards smaller numbers 
of stars (Eq. [5J). This can be shown as follows. Adopting 
for simplicity a mean stellar mass of m ~ 0.5 Mq we pre- 
dict (Eq. HJ that (dM/dt) dis ~ m(dN/dt) ~ m N/t dis ~ 
M 1_:E (ln0.02M/m) :E (Giersz & Heggie 1996) with x ~ 0.75 



6.3 The values of to c 



after core collapse 



The top right part of Fig. [5] shows the values of io° as a 
function of toM?, which is a proxy for the lifetime of the 
cluster. The pattern is the same as on the left side of this 
figure, i.e. before core collapse. The mean stellar mass at 
core collapse is shown in the middle right part of the figure. 
We find that we can approximate 



logmcc = +0.200 - 0.0984 x log toM ' 65 
logm cc = +0.075 - 0.0984 x log toM 080 



for Wo = 5 
for W (*2| 



The smaller mean mass of the Wo — 7 models at core col- 
lapse is due to the fact that core collapse occurs much earlier 
than for the less concentrated Wo — 5 models, so that the 
preferential loss of low mass stars had a smaller effect (see 
Fig.©. 

The lower right part of Fig. [5] shows the derived values 
of tr et . The mean values are log(t^ f cc ) = 0.864 ± 0.044 and 
0.796 + 0.017 for the models with Wo = 5 and 7 respectively. 
The scatter in t„f cc is larger than that of i„L f in Fig. \5\ 
because of the larger noise in the (dM /At) - vs - M relation 
due to the smaller numbers of stars (see Fig.rj}. We see that 
£^ f cc is slightly different for Wo — 5 and 7 models. This 
can be understood because inside a cluster the relaxation 
time increases with radius, so the outer cluster parts still 
keep some memory of the initial distribution by the time 
the center has gone into collapse. 
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7 CLUSTERS IN ELLIPTICAL ORBITS 

Models 21 to 25 are for clusters in elliptical orbits with a 
apogalactic distance of 8.5 kpc and various eccentricities. 
These can be compared with an otherwise similar model nr 

8 in a circular orbit. Fig. |8]shows the mass loss of models 8, 
21 and 23, with e = 0.0, 0.2 and 0.5 respectively. Since the 
total lifetime of the clusters scales approximately as (1 — e), 
the lifetimes of the clusters are very different (see Table Q]). 
The evolution can be described by the same three phases 
that we found for all other models: a stellar evolution domi- 
nated phase (A), a dissolution dominated phase before core 
collapse (B) and the phase after core collapse (C). The mass 
loss rates in all phases is variable with a periodicity of the or- 
bital period. The mass loss rate is highest at perigalacticon 
and the amplitude of the variations increases with increas- 
ing ellipticity. Especially after core collapse the amplitude 
increases drastically. This is due to the expansion of the 
outer layers of the cluster as a reaction to the core collapse. 
The stars in the outer layers are then more susceptible to 
the periodically changing tidal field. Only the models with 
e = and 0.2 show a jump in the mass loss rate at the time 
of core collapse. For clusters in more eccentric orbits no clear 
jump is observed at core collapse, but the mass loss rate does 
increase compared to the simple power-law extrapolation of 
phase B. 

The straight full lines in Fig. [8] show the values of 7 
defined by -(dM/dt)dis = Af 1 " 7 /t . Both the values of 7 
and of to describe the "time averaged" mass loss rate as a 
function of the mass. These values were derived from a study 
of the M(t) history for these models. The values of to, ^deiay, 
/™| x and ig c are listed in Tabled] 

Fig. [9] shows the values of to, m and t^ f of clusters 
in elliptical orbits (models 20 to 25) before and after core 
collapse respectively. We have added the data of model 8 
which has the same initial mass, Wo and the same i?Gai 
but a circular orbit. The values are compared with those of 
clusters with Wo = 5 at i?Gai = 8.5 kpc but with different 
masses (dashed lines), taken from Figs. [5] Notice that the 
models in eccentric orbits have very similar characteristics 
compared to those in circular orbits, if we correct for the 
effect of eccentricity by a factor (1 — e)" 1 in to and tg c . 
BM03 already concluded that the total life time of clusters 
is proportional to (1 — e). We find here that the "orbital 
averaged" mass loss rates before and after core collapse both 
scale with (1 — e) _1 . 

The main difference between the cluster models in cir- 
cular and elliptical orbits is in the mean stellar mass at core 
collapse (compare the middle panels of Figs. [S]and[5Jl. For 
clusters in circular orbits m cc increases towards shorter life- 
times, whereas m cc is about constant for clusters with the 
same apogalactic radius but different ellipticities. This is be- 
cause the BM03 models are Roche-lobe filling at perigalac- 
ticon so both t r h and t cc decrease steeply with increasing ec- 
centricity. The combination of a short lifetime (i.e. a higher 
maximum star mass at t cc ) and a smaller ratio of tcc/toM 7 
(i.e. fewer low mass stars are lost after mass segregation) 
results in m C c being almost independent of e. This is illus- 




4.2 4.0 3.8 3.6 3.4 3.2 3.0 

log(M) [MJ 



Figure 8. The mass loss rates of clusters with elliptical orbits 
with eccentricities of (upper), 0.20 (middle) and 0.5 (bottom). 
A = dominated by stellar evolution, B = dominated by dissolu- 
tion, C = dominated by dissolution after core collapse. The line 
separating phases B and C indicates the core collapse time. The 
figure shows (dAf/dt) versus M in logarithmic units. The full up- 
per line is the total mass loss rate; the dashed line is the mass loss 
by stellar evolution; the full smooth lines are the mass loss by dy- 
namical effects, assumed to scale as a power law of M , before and 
after core collapse. The models are defined by a vector contain- 
ing: model nr, nr of stars, total lifetime (in Gyrs), concentration 
parameter Wo, Rg (in kp c )i and orbit. 



trated in Fig. 1101 which shows the variation of m with /x and 
the mean mass at core collapse, which is almost constant. 

The core collapse time of the clusters in elliptical orbits, 
shown in Fig. [7] is slightly longer than predicted by Eq. [20] 
by about 20 percent. This is because the values of t r h, listed 
in Table [T] refer to perigalacticon, which is smaller than the 
orbital-averaged values of t r h. 
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Figure 9. The parameters of the models of clusters in elliptical orbits, nrs 20 to 25, combined with model 8 of the same mass but with 
e = before (left) and after (right) core collapse. In all panels the right data point is for e = and the left one is for e = 0.8. Top: The 
relation between toM^ and to before (left) or tg c (right) after core collapse. The dashed lines show the relation for clusters in circular 
orbits at Rq^i = 8.5 kpc (models 6 to 10) for comparison (Fig. |5J. Middle: The mean stellar mass before (left) and at core collapse 
(right), again compared to the that of models 6 to 10. Bottom: The resulting values of tfv, (left) and t^ f cc (right) as a function of tgM? . 



The mean values of tiL* 



13.3 Myr and ri. cc = 7.2 Myr are the same as for models 6 to 10 with circular orbits. 




-0.4 -0.6 -0.8 -1.0 

log ij. = log M(t)/M, 

Figure 10. The variation of the mean stellar mass as function 
of ii for clusters in elliptical orbits. Lowest curve: e = 0.2, upper 
curve: e = 0.8. The wiggles are the result of the variation in the 
mass of the cluster within the variable tidal radius. Diamonds 
indicate the minimum values and squares indicate the moment of 
core collapse. Notice that the mean mass at core collapse is about 
the same for all models. 



INITIALLY ROCHE-LOBE UNDERFILLING 
CLUSTER MODELS 



8.1 The parameters of the initially Roche-lobe 
underfilling models 

Because we are also interested in the mass history of initially 
highly concentrated clusters, i.e. with a half-mass radius 
much smaller than their tidal radius, we have supplemented 
the set of Roche-lobe filling cluster models with TV-body 
simulations of a series of clusters that start severely Roche- 
lobe underfilling. These are for clusters with an initial mass 
in the range of 10400 to 84000 M@, in circular orbits at 
-RGai = 8.5 kpc and with an initial density distribution de- 
scribed by a King parameter of Wo — 5. The metallicity is 
Z = 0.001. These cluster models have an initial half mass 
radius of r^ = 1, 2 or 4 pc. The stellar IMF of these clusters 
is different from those of the Roche-lobe filling model. They 
have a Kroupa mass function in the range of 0.10 to 100 
Mq, with an initial mean stellar mass of 0.623 Mq. In these 
models 10% of the formed neutron stars and black holes are 
retained in the cluster. The models were calculated for this 
study. The parameters of the models are listed in Table [2] 

For the study of the mass loss by dissolution we define 
an "underfilling factor" 5, which is defined as the ratio of 
the initial half- mass radius of the cluster model, rh, and the 
half-mass radius r{j f of a Roche-lobe filling cluster with the 
same King parameter Wo- 
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3w Q 



£h_ _ (VhAt)wo 



''.7 



(23) 



where r j = r\ is the Jacoby radius, i.e the tidal radius of a 
Roche-lobe filling cluster, and r t ° is the end of the density 
profile of a cluster with fixed values of r^ and Wo- Note that 
5=1 for Roche-lobe filling clusters and < 1 for Roche-lobe 
underfilling clusters. In this expression (fh/rt)s = 0.187 and 
(rh/rt)r = 0.116 are the ratios for a Roche-lobe filling clus- 
ter with an initial density distribution of Wo = 5 and 7 
respectively. (The same definition of 5 was used by Gieles 
an Baumgardt (2008) in their theoretical study of the mass 
loss of Roche-lobe underfilling clusters.) The values of 5s are 
listed in Tabled We will show below that for the description 
of the dissolution of the Roche-lobe underfilling models the 
parameter 3V = 1.6125s is more important than 5s. 

8.2 Dissolution before core collapse 

The mass loss of the Roche-lobe underfilling models can be 
described in the same way as for the Roche-lobe filling mod- 
els. The resulting fitting parameters are listed in the right 
hand block of Tabled 

The first phase is dominated by stellar evolution. How- 
ever in this case there is no evolution-induced mass loss. 
This is because the clusters are initially well within their 
tidal limit. So the mass loss by stellar evolution does pro- 
duce an expansion of the radius, but this expansion does 
not immediately reach the tidal radius. This is reflected in 
the values of f^ x = 0.0 for all models, except uf3 and uf6. 
These two are the least Roche-lobe underfilling models with 
5s — 0.66 and 0.53 respectively. The value of /™ can then 
be described in a similar way as for Roche-lobe filling clus- 
ters, Eg. 1161 but with a correction term that depends on the 
underfilling factor, 



/Sd x = -0.86 + 0.40 x log(ioMr) + 2.75 x log & 



(24) 



with a maximum of 1 and a minimum of and 7 = 0.80 (see 
below) . 

The dissolution before core collapse can be expressed by 
a power law approximation of (dM/dt)^° s rm versus M, with 
(dM/dt)3° rm defined by Eq.[T7| The lower left panel of Fig. 
[Jshows the relation between log((dM/d£)d° rm ) and log (M). 
Because the core collapse time of the Roche-lobe underfill- 
ing models is short, each model contributes only to a small 
part of the mass range. We found that 1 — 7„f = 0.20 ± 0.05 
gives a good fit to these plots, so we adopted 7 u f = 0.80. We 
note that this value is the same as for the Roche-lobe filling 
clusters with Wo = 7, whereas the initially Roche-lobe un- 
derfilling clusters have an initial concentration of Wo = 5. 
However, a study of the expansion of severely Roche-lobe 
underfilling clusters by means of iV-body- simulations has 
shown that the initial expansion due to mass loss redis- 
tributes the density close to that of a Wo = 7 model, in 
agreement with our derived value of 7 u f = 0.80. 

The tendency of clusters to evolve to a Wo ~ 7 model 
was noticed by Portegies Zwart et al. (1998). In initially 
strongly concentrated clusters, Wo > 7, dynamical friction 
quickly drives the massive stars to the center where they 



will lose mass due to stellar evolution. This results in an 
expansion of the core and a less steep density profile, so Wo 
decreases. On the other hand, in clusters with a less steep 
initial concentration, Wo ^ 5, dynamical friction is less ef- 
ficient and the massive stars lose mass by stellar evolution 
before they reach the center. So the cluster expands more 
homogeneously due to evolutionary mass loss. After the mas- 
sive stars have undergone stellar evolution, dynamical effects 
take over and the core shrinks due to the approaching core 
collapse, so the density distribution becomes more concen- 
trated and Wo increases. 

Once the cluster has expanded to about the tidal radius, 
the dissolution is very similar to that of a Roche-lobe filling 
cluster with Wo = 7. This is reflected in the values of £q 
which range from 5.4 to 7.5 Myr, whereas the values of to 
for the comparable Roche-lobe filling models, nrs 16 to 19, 
range from 6.0 to 6.5 Myr. Part of the difference is due to 
differences in the mean stellar mass, because the two sets of 
models have different IMFs. 

Dissolution of the Roche-lobe underfilling models needs 
more time to get started than Roche-lobe filling models, be- 
cause the clusters first have to expand to the tidal limit. This 
is a slow process that occurs on the relaxation time scale. 
This is shown in Fig. [TT] which shows the ratio tdeiay/trh as a 
function of the Roche-lobe underfilling factor Js, for differ- 
ent values of r^. The values of tdeiay/trh range from about 
to 4 for the models considered here. The figure shows that 
we can approximate 



£dclay 



4.31 10' 



(8s 



-1.989 . 1.605 
trh 



(25) 



with all ages in Myr. This equation is valid for —1.0 < 
Iog^s) < —0.20. For Roche-lobe filling clusters with 
log (5s) > —0.20 the delay time does not scale with t r h 
but with the crossing time at the tidal radius and ideiay — 
3.0 tcr(rt) (Sect. [4j| . So, for strongly Roche-lobe underfilling 
models the delay time scales with the initial value of t r h, but 
if the underfilling factor approaches 5=1 the delay time is 
much shorter and scales with t CT (rt)- 

The dissolution before core collapse can now be ex- 
pressed by Eq. Q20 with /delay given by Eq. [T3] and tdoiay 
approximated by Eq. [25] 

Fig. [12] (left) shows the values of if! 1 t ne mean mass 
before core collapse m and the resulting values of t^f.uf as 
a function of toM^ 1 . We used different symbols for different 
initial half mass radii. 

The mean stellar mass before core collapse depends 
both on toM^ and on J and not only on toM^ as is the case 
for the Roche-lobe filling clusters. This is because the core 
collapse time t cc depends on 5 and so the amount of mass 
that the cluster has lost before core collapse also depends on 
both toM^ and t r h. We find that we can express the mean 
stellar mass before core collapse, i.e. between t(n = 0.25) 
and t cc , as 



logm = 0.139 -0.09841og(i M: y ) + 0.101 log(3V 



(26) 



We have expressed m in terms of $7 (with ^7 = 1.6125s, 
Eq. 123ft instead of the initial value of 5s because the initial 
expansion of the clusters redistributes the density to about a 
Wo = 7 model (see above). The relation has about the same 
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Table 2. The N-body models of initially Roche-lobe underfilling clusters. The left block of contains the model parameters, the middle 
section the cluster time scales and the right block contains our fit parameters. 



Input parameters 
# Mass nr 
Mq stars 


Wo 


kpc 


pc pc 


55 


Timescales 
Gyr Gyr Gyr 


t„(r h ) t 
Myr 


= r (n) 
Myr 


Output parameters 

7 t t M:< 

Myr Gyr 


Myr 


f& 


i 


Myr 


Jcc 


ufl 10831 


16k 


5 


8.5 


32.6 1.00 


0.164 


7.22 3.00: 0.045 


0.57 


75.5 


0.80 5.1 8.6 


50 


0.0 


10.5 


1.09 


uf2 10426 


16k 


5 


8.5 


32.2 2.00 


0.332 


7.59 3.50: 0.127 


1.65 


75.5 


0.80 6.2 10.2 


100 


c 


.0 


10.4 


1.34 


uf3 10589 


16k 


5 


8.5 


32.4 4.00 


0.660 


5.89 5.30: 0.361 


4.64 


75.5 


0.80 5.0 8.3 


100 


( 


.2 


8.4: 


1.15 


uf4 21193 


32k 


5 


8.5 


40.8 1.00 


0.131 


11.42 3.36 0.058 


0.11 


75.5 


0.80 5.0 14.5 


200 


c 


.0 


11.1 


1.16 


uf5 21095 


32k 


5 


8.5 


40.7 2.00 


0.263 


13.40 7.00 0.165 


1.16 


75.5 


0.80 6.0 17.3 


350 


c 


.0 


12.8 


1.11 


uf6 20973 


32k 


5 


8.5 


40.7 4.00 


0.526 


12.75 9.30 0.466 


3.29 


75.5 


0.80 6.0 17.2 


200 


c 


.1 


10.fi 


1.21 


uf7 41465 


64k 


5 


8.5 


51.0 1.00 


0.105 


17.79 3.73 0.076 
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Figure 12. The parameters of the models of initially Roche-lobe underfilling clusters, nrs ufl to ufll before (left) and after core collapse 
(right). Top: The relation between toM^ and tj} (left) or £q C ' u (right), i.e. after core collapse. Middle: The mean stellar mass before 
(left) and at core collapse (right). The dashed lines indicate mean relations of the different values of r n . Bottom: The resulting values of 
*ref uf (teft) an d Ket cc uf ( r 'S nt ) as a function of toM7. The mean relations are given in Sect. 18. 2l and 18.31 



dependence on toM^ as in the case of Roche-lobe filling 
clusters (Eq. [i"8)) . 

The lower left part of Fig. [12] shows the resulting values 
of i^ f uf . For Roche-lobe filling clusters we found one value 
of i^ f = 3.5 Myr (logi£ f = 0.544) for all models of W = 
7 (see Fig. [5|. The Roche- lobe underfilling models show a 
large range in t^ ci uf , indicating that the underfilling factor 
plays is a role for small values of $7. We found that we can 
approximate 



log(i 



N \ 
rcf ,uf / 



= 0.544 if log(g T ) > -0.50 

0.544 + 0.65 x {log(5 7 ) + 0.50} 
if log($ 7 ) < -0.50 



(27) 



This shows that there is a smooth transition between the 
dissolution parameter for Roche-lobe filling and underfilling 
clusters. The fits are shown as three partially overlapping 
dashed lines for models of 7- n =l, 2 and 4 pc in the lower 
part of Fig. [H 



8.3 Dissolution after core collapse 

The time of core collapse is expected to scale with the half 
mass relaxation time. Fig. [13] shows t cc as a function of i rn 
and the initial Roche-lobe underfilling factor r\/rt. In Sect. 
I6.1l we found that t cc oc £ r h°' 872 for Roche-lobe filling clus- 
ters. We find that for Roche-lobe underfilling clusters the 
dependence of t cc on i rn can be described by the same power 
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initially Roche-lobe underfulling clusters. The dashed lines show 
the fit of Equ. [25] 



Figure 13. The core collapse time of the Roche-lobe underfilling 
models as a function of the initial half mass relaxation time and 
the underfilling factor 5s- 



law dependence, but there is an additional dependence on 
the underfilling factor 

log(tcc) ^ 1.505 + 0.872 log(t rh ) - 0.513 log (5s) (28) 

(In this case ^5 is the crucial parameter, rather than 3V, be- 
cause the initial relaxation time t r h is defined for the initial 
density distribution with Wo — 5.) The smaller the under- 
filling factor, the larger t cc for a given value of t r h. This is 
because clusters that are initially strongly Roche-lobe un- 
derfilling expand more strongly and so t r h(t) increases more 
strongly with time, which results in a larger ratio of t cc /t T h- 
In the limit of $5 = 1, i.e. for Roche-lobe filling clusters Eq. 
[281 predicts that log(t cc ) ~ 1.50 + 0.872 log(i rh ). This can 
be compared with the value for Roche-lobe filling clusters 
of log(i cc ) ~ 1.23 + 0.872 log(trh) (Eq. H]}. So Roche-lobe 
underfilling clusters need about twice as many "initial" re- 
laxation times to go into core collapse as Roche-lobe filling 
clusters. This is because of their stronger expansion and the 
fact that they do not push stars over the tidal boundary 
while evolving towards core collapse. 

A study of the relation between (dM/d£)5° rm and M(t) 
of the various Roche-lobe underfilling models shows that the 
dissolution after core collapse can be described by the same 
values of 7 CC as for the Roche-lobe filling cluster models, i.e 
7cc = 0.70 if M(t) > 10 3 M Q and 7cc2 = 0.40 if M(t) < 10 3 
Mq (see lower right panel of Fig. [4j . This is due to the fact 
that core collapse results in a redistribution of the density in 
the cluster and erases the memory of the pre-collapse phase. 

The values of ig c of the Roche-lobe underfilling mod- 
els are listed in Table [2] They are plotted versus toM? in 
the right hand panel of Fig. [T5] together with the mean 
stellar mass at core collapse and the resulting values of 
t^ f cc . The dissolution parameter after core collapse when 
M(t) < 1O 3 M is £g c2 = tff x io 3 <7cc-7cc2) = 10 o.90 f cc 

The mean mass at core collapse can be approximated 
by 

log(m cc ) = 0.178 - 0.0984 log(t M^) + 0.207 log(3V) (29) 



The relation has the same slope as that of the Roche-lobe 
filling clusters after core collapse (Eq. I22[l but the constant 
0.178 for clusters with $7 = 1 is different from the 0.075 of 
Roche-lobe filling clusters with Wo — 7, which indicates a 
higher mean mass at core collapse. This is mainly due to 
differences in the IMF of the two sets of models. 

Fig.[T2p shows that t* t yCC = t„ c xm cc J with 7 = 0.70 for 
these models is about constant but with a significant scatter. 
This shows that there is a small residual effect of the initial 
underfilling factor in the dissolution after core collapse. We 
found that we can approximate 



iog(t£rf )CC ) 



0.875 + 0.127 x log(5 7 



(30) 



The constant 0.875 is slightly larger than the value of 0.796 
for Roche-lobe filling models of Wo — 7. 

Fig. [14] shows the history of the mean stellar mass in 
the Roche-lobe underfilling cluster models as a function of 
fj,. The trends are approximately the same as in Fig. [6] i.e. 
an initial decrease due to the loss of massive stars by stellar 
evolution followed by an increase of m after mass segregation 
has been established and low mass stars are lost preferen- 
tially. The difference is due to the different initial stellar 
IMF, which have an initial m = 0.623 M instead of 0.547 
for the Roche-lobe filling models. Moreover, in the Roche- 
lobe underfilling models 90% of the black holes and neutron 
stars is ejected whereas they were all retained in the Roche- 
lobe filling models. Both effects influence the evolution of rn 
(Kruijssen 2009). The decrease of m is stronger than for the 
Roche-lobe filling models because the IMF of the Roche-lobe 
underfilling models reaches up to 100 Mq , whereas the IMF 
of the other models reaches to 15 Mq . We see that m of the 
Roche-lobe underfilling models reaches its minimum at the 
same value of t ~ 0.15itot as the other models. This shows 
that full mass segregation is reached at the same fraction of 
the total lifetime, independent of the initial radius. 
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9 THE RELATION BETWEEN THE CLUSTER 
LIFETIME AND t M? 

The lifetime of clusters depends on time scales for mass loss 
by stellar evolution, and the dissolution constants to and tg c 
before and after core collapse. Stellar evolution dominates 
the mass loss only during the first part of the cluster life- 
time and typically removes about 20 to 40% of the cluster 
mass, after which dissolution takes over. Since the dissolu- 
tion timescales before and after core collapse both depend 
on the strength of the tidal field in which the cluster moves, 
we may expect that the cluster lifetime depends largely on 
toMC- 

This is confirmed in Fig. 1151 which shows a very tight 
relations between to x MC and t\%. For clusters with an 
initial density concentration described by Wo = 5 and Wo = 
7 in circular and elliptical orbits we find 



log(ti%) = 0.518 + 0.864 x log(t Mi°' 6j ) if W = 5 

= 0.797 + 0.778 x log(t Mi°' 80 ) if W = 7(31) 

The relation between toM? (with 7 = 0.80) and t\% for 
the Roche-lobe underfilling clusters is indistinguishable from 
that of the Roche-lobe filling clusters of Wo = 7. The tight 
correlations show that toM^ can be used as accurate indi- 
cator of the lifetime of a cluster. 

Equation 1311 might suggest that the lifetime of a cluster 
is proportional to A/i ' 56 for W = 5 clusters and Mi ' 62 for 
Wo = 7 cluster. However, we remind that to is proportional 
to m-~< (Eq.HD and m oc (toMj 7 )" - 121 and (ioM^)" 0094 
(Eq. I18[) for Wo = 5 and 7 respectively. This implies that 
t x % oc Mi 061 and Mi 067 for W - 5 and 7, in agreement 
with the values of the indices 0.62 and 0.67 derived by BM03. 



10 THE PREDICTED MASS HISTORY OF 
STAR CLUSTERS 

In the previous sections we have described the interplay 
between the different mass loss processes of star clusters. 
Based on these results we derived a recipe for calculating 
the mass evolution of star clusters in different environments. 
The recipe is described in Appendix A. 

10.1 The contribution of different effects to the 
mass loss 

The mass loss from star clusters is due to several effects: 
stellar evolution, evolution-induced loss of stars, and dis- 
solution (relaxation-driven mass loss) before and after core 
collapse. Fig. [16] shows the contributions of these different 
effects for two characteristic models, # 15 which has a life- 
time of ti% = 1.3 Gyr and # 2 with t 1% = 26.9 Gyr. The 
two models show that clusters with a long lifetime (> 20 
Gyr) lose about 35% of their mass by stellar evolution, 15% 
by induced mass loss and the remaining 50% by dissolution. 
Clusters with a short lifetime (< 5 Gyr) lose more than 60% 
by dissolution, less than about 30% by stellar evolution, and 
less than 10% by induced mass loss. This is because the short 
lifetime is the "result" of a strong mass loss by dissolution, 
which does not leave much time for the cluster to lose a large 
fraction of its mass by stellar evolution. 

10.2 Predicted M(t) of the total mass of clusters 

We have calculated the M(t) history of all cluster models 
listed in Table Q] with the recipe described in Appendix A. 
A subset of the results is shown in Fig. 1171 The sample 
shown contains models with Wo = 5 with a large (64k or 
128k) and small (8k) number of stars at respectively Rg^i = 
15, 8.5 and 2.63 kpc; two models with W = 7 (128k and 
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8k); two models in elliptical orbits (e = 0.2 and 0.5); and 
five Roche-lobe underfilling models with different numbers 
of stars and initial half mass radii. The agreement is good for 
all models of clusters in circular and elliptical orbits, with 
initial concentrations Wo — 5 and 7, and for the Roche-lobe 
underfilling clusters, including the ones not shown here. For 
cluster models in the original BM03 sample that are not 
discussed in this paper, the agreement is equally good. 

The different models have different shapes of M(t)/Mi 
versus t/ti%. Clusters with a long lifetime (ti% > 20 Gyr, 
high Mi) show a strong drop in mass during the first 5% 
of their life, due to stellar evolution and induced mass loss, 
followed by a more gentle decrease. Clusters with a short 
lifetime (t 1 % < 10 Gyr, low Mi) show a more gradual con- 
cave shape. All models show a bump in the M(i)-plot near 
the core collapse time: the mass loss rate is about twice as 
high after core collapse than before. 

The shapes of the M (t) relations are all convex with 
various degrees of curvature. Only those cluster models for 
which core collapse occurs about halfway through their life- 
time show a more or less linear mass history (e.g. models 
16, uf09 and uflO). 



11 DISCUSSION 

We have shown how the different mass loss effects of star 
clusters interact in the determination of their mass history. 
We have also derived a recipe for calculating the mass loss 
history of star clusters in different environments and with 
different metallicities and stellar initial mass functions. This 
study is based on the iV-body simulations by BM03, supple- 
mented with newer iV-body simulations of Roche-lobe un- 
derfilling clusters, so the results are dependent on the char- 
acteristics of these models. The N-boAy models that we used 
are relatively simple: the clusters start in virial equiliblium, 
without primorial mass segregation, without primordial bi- 
naries and with stars in isotropic orbits. The study of these 
simple models, which are valid after the gas expulsion phase, 



are a first step to understand the complicated interplay be- 
tween the various dynamical effects in clusters. The mod- 
els can be refined later when observational evidence shows 
which assumptions have to be improved. 

We discuss the major assumptions of the models, how 
they may have influenced our results, and how they can be 
taken into account in the recipe for computing M(t). 

(a) The BM03 models are Roche-lobe filling. This is 
a good assumption for open clusters and globular clusters 
which are close to the galactic centre or have very large 
half-mass radii. However, the majority of globular clusters 
probably formed with half-mass radii around 1 pc and there- 
fore started strongly Roche- lobe underfilling (Baumgardt, 
Kroupa & Parmentier 2008b). Some of these are still under- 
filling at present (Baumgardt et al. 2010) 

(b) The models that we used do not have initial mass 
segregation. The question of the initial mass segregation is 
still open. Baumgardt et al. (2008a) found that the present 
overall mass function of most globular clusters can be ex- 
plained without invoking initial mass segregation, but some 
clusters require initial mass segregation to explain their 
present mass function. Observations of young clusters, age 
<50 Myr, show evidence for mass segregation, e.g. Brandl 
et al. (1996) for R136; Hillenbrand and Hartmann (1998) for 
the Orion Nebula Cluster; McGrady et al. (2005) for M82- 
F. However, de Grijs et al. (2002) argued that this does 
not necesaarily imply initial mass segregation because the 
timescale for the dynamical segregation of high mass stars in 
young massive clusters may be very short. Dynamical mass 
segregation has been taken into account in the models we 
used in this study. 

(c) The models have no primordial binaries, but only 
dynamically formed binaries. It is expected that real clusters 
contain a large fraction of initial binaries (Elson et al. 1998, 
Hut et al. 1992, Hu et al. 2006, Sommariva et al. 2009). 
However only hard binaries influence the cluster dynamics 
as they can heat the cluster and prevent core collapse (Hut 
et al. 1992). Kiipper et al (2008) have shown that the escape 
rate is hardly affected by binaries. 

(d) Stellar remnants are initially retained in the BM03 
cluster models, i.e. they are not ejected by a kick-velocity. 
The new Roche-lobe underfilling models have a 10% reten- 
tion factor of black holes and neutron stars. The retention of 
the remnants implies that the model clusters at older ages 
contain a large fraction of neutron stars which are more 
massive than the average stellar mass. Clusters with a large 
fraction of massive remnants will dissolve faster due to the 
higher average stellar mass. Also the depletion rate of low 
mass stars could be different (Kruijssen 2009). 

(e) The effects of bulge shocks are included in the mod- 
els of clusters in eccentric orbits. However, disk-shocking is 
not included in our models. Vespirini and Heggie (1997) have 
studied the effect of disk-shocking. Based on their results 
(see their Fig. 21) we conclude that the effect of disk-shocks 
on decreasing the lifetime of clusters is small, especially for 
clusters beyond the solar circle and for massive clusters. Our 
models also do not include the effects of shocks by spiral 
density waves or passing GMCs. The latter effect is thought 
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Figure 17. The M(t) history of a representative set of cluster models: 6 models in circular orbits with Wo = 5, two with Wo = 7, two in 
elliptical orbits, and 5 Roche-lobe underfilling models. The models of the Roche-lobe filling clusters are specified by a vector containing: 
model nr, number of stars, Wo, i?Gal m kpc and orbit. The models of the Roche-lobe underfilling clusters are specified by: model nr, 
number of stars, Wo, Rgs.1 m kpc and r^ in pc. Full lines: derived from TV-body simulations of BM03 or from the sample of Roche- lobe 
underfilling models presented here. Dotted line: predicted with the parameters listed in Tab.[T] Dashed lines: predicted by the method 
described in Appendix A. The upper dash-dotted line shows the fraction of the mass that is lost directly by stellar evolution. The vertical 
tickmarks indicate the time of core collapse from the models (full) and calculated with Eqs, 1201 and 1281 (dashed). 



to be the main destruction mechanism for clusters in the 
Galactic plane (e.g. Lamers & Gieles 2006) and probably 
also in GMC-rich interacting galaxies (Gieles et al. 2008). 
These effects will increase the dissolution compared to that 
of the BM03 models. However they can easily be accounted 
for in the recipe that we derived for calculating M(t), by 
simply adopting a smaller value of the dissolution parame- 
ter to- This correction is justified because shocks will remove 
stars from the outer regions of the cluster in approximately 
the same way as the tidal field. 



(f) The models are calculated for a given stellar initial 
mass function, (a Kroupa mass function of 0.15 < m < 15 
Mq for the BM03 models and a Kroupa mass function of 
0.10 < m < lOOM© for the new Roche-lobe underfilling 
models) and for a given metallicity of Z — 0.001. The mass 
loss by stellar evolution depends on these assumptions. How- 
ever, the recipe that we derived allows the choice of different 
metallicities and different IMFs, by applying the approxi- 
mate formulae that describe the mass loss by stellar evolu- 
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tion and the formation of remnants for a grid of metallicities 
listed in Appendix B. 

(g) We assumed that the cluster move in a spherical 
logarithmic potential with a constant rotation speed. This 
implies that we may have underestimated the effect of disk- 
shocking, which is important for clusters in disk galaxies. A 
study of the effects of a non-spherical halo and the resulting 
non-circular orbits with disk-shocking has to be postponed 
to future studies. (The use of GPUs for the computations 
of cluster dynamics will allow a significant expansion of the 
parameter space of cluster models.) 



12 SUMMARY AND CONCLUSIONS 

Based on iV-body simulations by BM03 of the evolution of 
Roche-lobe filling star clusters of different initial concentra- 
tions and in different orbits in the Galaxy, and on a new 
sample of Roche-lobe underfilling clusters in circular orbits, 
we have studied the interplay between the different mass loss 
effects: mass loss by stellar evolution, loss of stars induced 
by stellar evolution, and dynamical mass loss (referred to as 
"dissolution") before and after core collapse. 

At young ages stellar evolution is the dominant effect. 
The fast (adiabatic) evolutionary mass loss results in a si- 
multaneous expansion of the cluster and a shrinking of its 
tidal (Jacobi) radius. So the outer cluster layers become un- 
bound. This evolution-induced mass loss contributes to the 
overall mass loss if the cluster is deeply emersed in the tidal 
field, i.e. if the cluster is initially filling its tidal radius. The 
evolution-induced mass loss rate is proportional to the mass 
loss rate by stellar evolution but it is smaller for clusters 
with a mass loss rate by dissolution larger than the mass 
loss rate by stellar evolution. This is for instance the case 
for low mass clusters or for clusters in orbits close to the 
Galactic center. The iV-body models show that the induced 
mass loss does not start immediately, but that it needs time 
to build up. This build-up can be described by an expo- 
nential function (Eqs. [13] and I14|l with a delay time scale 
of a few, typically 3, times the crossing time at the tidal 
radius. For Roche-lobe underfilling clusters the delay time 
scale is much longer, of the order of a few half mass relax- 
ation times, because the cluster first has to expand to the 
tidal radius. The actual value depends on the initial Roche- 
lobe underfilling factor. As the evolution-induced mass loss 
rate needs time to get going, the total amount of evolution- 
induced mass loss is considerably smaller, typically 10 to 
50% of the total amount of mass lost by stellar evolution 
(see Fig. Dpi , 

The mass loss of the cluster models by dissolution needs 
time to build-up, just like the evolution-induced mass loss. 
However, this is a consequence of the initial conditions of 
the cluster model and the start of the dissolution might be 
very different in real clusters (Sect. 5.3). 

We have shown from both theory and the model simula- 
tions that the dissolution rate depends on the environment 
of the clusters and can be described accurately by a formula 
of the type (dM/di) d i s = -M (i) 1_7 /t , with M and t in 
units of solar mass and Myr. The value of the dissolution 



parameter to depends on the environment, e.g. the Galactic 
potential, the orbit and shocks by spiral arms and passing 
GMCs. 

We have derived expressions for estimating to for clus- 
ters in galaxies where tidal evaporation is the main dissolu- 
tion effect. This value depends on the Galactic potential (i.e. 
the galactic rotation velocity), the orbit of the cluster, the 
initial concentration characterized by Wo, and on the evo- 
lution of the mean stellar mass. We have derived an expres- 
sion for the mean stellar mass during the pre-core collapse 
phase and the post-core collapse phase for various initial 
mass functions and metallicities. For clusters in an environ- 
ment where shock-heating by encounters with spiral arms or 
GMCs are important, the value of to can be estimated using 
the descriptions by Gieles et al. (2006, 2007). 

The slope of the stellar mass function depends mainly 
on the remaining mass fraction /j, = M(t)/M{ of the clusters 
and hardly on the initial parameters, such as mass, mass 
function, concentration factor and strength of the tidal field 
(e.g. Vesperini & Heggie, 1997; BM03; Trenti et al. 2010). 
This effect can also be seen in our results in the evolution 
of the mean stellar mass of a cluster. The data in Figs. [6] 
1101 and [14] show a very similar evolution of m as function of 
H in almost all models. The main difference is that clusters 
with short lifetimes have an offset of m to higher values. 
This is because the mean stellar mass not only depends on 
the slope of the mass function, but also on the mass of the 
most massive stars that have survived stellar evolution. This 
upper mass depends on the age of the cluster and not on its 
mass fraction. 

The details of the iV-body simulations have shown that 
7 = 0.65 for clusters with an initial density distribution 
of a King-profile with Wo = 5 and 7 = 0.80 if Wo = 7. 
The difference in 7 is due to the fact that the dissolution 
timescale depends on both the half-mass relaxation time and 
the crossing time. Initially Roche-lobe underfilling clusters 
quickly expand due to mass loss by stellar evolution and 
reach a density distribution of approximately Wo = 7, so 
their dissolution is also described by 7 = 0.80. These values 
of 7 apply to the pre-core collapse phasqj. 

We note that our Roche-lobe underfilling models have 
half-mass relaxation times between 40 and 800 Myr. The 
central relaxation times are about 10 times shorter, i.e. 4 
to 80 Myr, but still longer than the evolution time of the 
most massive stars. If the initial radius is smaller than those 
of our models, e.g. ^ 0.5 pc, the core relaxation time may 
be shorter than the evolution time and the cluster concen- 
tration might decrease rather than increase due to stellar 
evolution. 

The iV-body-simulations showed that cluster dissolu- 
tion does not start right away but that it also needs time 
to get going. We find that this build-up can be described by 
the same exponential function with the same time scale as 
the evolution-induced mass loss (Eq. ll4|) . 

4 We point out that this formula describes the time-dependent 
mass loss rate per cluster. It is different from the formula that 
was derived by BM03 (their Eq. 7) to describe the dependence 
between the total lifetime of a cluster and its initial mass. 
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The core collapse time t cc of the models can be ex- 
pressed in terms of the initial half mass relaxation time £ r h 
by a simple relation that depends on the underfilling factor 
5 (Eq.HSJl. For a Roche-lobe filling cluster of M = 10 5 M© 
in a circular orbit i cc ~ 5 t Th (t = 0). 

The mass loss rate by dissolution increases at core col- 
lapse by about a factor 2 depending on the model, and has 
a different mass dependence after core collapse than before 
with (dM/dt)dis = -M{t) x ^ cc /tg c with 7cc = 0.70 for all 
models. This is independent of the initial density distribu- 
tion because this is erased by the core collapse. When the 
mass of the cluster decreases to M(£) ^ 10 3 M© the mass 
loss dependence changes to 7 CC = 0.40. This is due to the 
variation of the Coulomb logarithm in the dependence of re- 
laxation time on the number of stars ./V in the cluster. We 
derived an expression for tg c , i.e. after core collapse (Sect. 

E3J. 

We have derived simple expressions for the parame- 
ters that describe the evolution-induced mass loss and the 
dissolution, in terms of the initial cluster parameters (Mi, 
Wo) and the orbit (i?Gai and eccentricity). We also de- 
rived parameters that describe the mass loss by stellar 
evolution for different stellar IMFs and metallicities. With 
these parameters we can describe the different mass loss ef- 
fects throughout the lifetime of a cluster. By integrating 
(dM/dt)tat = (dM/di)ev + (dM/dt)%i + (dM/d£) dis , start- 
ing from the initial mass Mi, we can calculate the mass loss 
histories of clusters. For this purpose we describe a simple 
recipe for calculating M(t) in Appendix A, that provides a 
summary of the equations. The resulting mass histories are 
compared with those derived from the iV-body simulations. 
Some of the characteristic results are shown in Fig. 1171 The 
agreement is very good, within a few percent of the initial 
mass. The agreement is equally good for the cluster models 
of BM03 that were not used in this paper. 

The method described here provides a description of 
the variation of the total mass, i.e. stars and remnants, of 
a cluster with age. To derive the luminous mass, one has to 
correct the total mass for the contribution by remnants. In 
the calculations of BM03 the newly formed remnants were 
retained in the cluster (no kick velocity was assumed). In 
the method that we present here the kick fractions of black 
holes, neutron stars and white dwarfs can be specified as 
free parameters. In later phases part of the remnants can 
be lost by dissolution. At late ages remnant neutron stars 
and black holes are the most massive objects in the cluster. 
They will sink to the center and are not likely to be lost by 
dynamical effects. 

The results of this paper and the methods can be used 
to predict the mass histories of star clusters with different 
stellar IMFs and different metallicities in different environ- 
ments. This can then be used to predict the evolution of the 
mass function of cluster systems. 



ACKNOWLEDGMENTS 

We thank Onno Pols for providing us with updated evo- 
lutionary calculations and Diederik Kruijssen for his com- 



ments on the manuscript. Simon Portegies Zwart has given 
important advice. HJGLM and HB thank ESO for Visiting 
Scientist Fellowship during a few moths in 2008 and 2009 
in Garching and Santiago, when this study was performed. 
This research was supported by the DFG cluster of excel- 
lence Origin and Structure of the Universe (www. universe- 
cluster. de). 



REFERENCES 

Aarseth S.J., 1999, PASP, 111, 1333 
Aarseth S.J., Lecar M., 1975, ARAA, 13, 1 
Ambartsumian V.A. 1938, Ann. Leningrad State Univ., No 

22, Vol 4, p. 19 
Anders P., Fritze-v. Alvensleben U., 2003, A&A 401, 1063 
Baumgardt H., 2001, MNRAS, 325, 132 
Baumgardt H., de Marchi G., Kroupa P., 2008a, ApJ, 685, 

247 
Baumgardt H., Kroupa P., Parmentier G., 2008b, MNRAS, 

384, 1231 
Baumgardt H., Makino J., 2003, MNRAS, 340, 227 (BM03) 
Baumgardt H., Parmentier G., Gieles M., Vespirini E., 

2010, MNRAS, 401, 1832 
Boutloukos S.G., Lamers, H.J.G.L.M., 2003, MNRAS, 338, 

717 
Brandl et al. 1996, ApJ, 466, 254 
Bruzual G., Chariot, S., 2003, MNRAS, 344, 1000 
Chandrasekhar S., 1943, ApJ, 98, 54 
Chernoff P., Weinberg M., 1990, ApJ, 351, 121 
De Grijs R., Johnson R.A., Gilmore G.F., Frayn CM., 

2002, MNRAS, 331, 228 
Elson R.A.W., Sigurdsson S., Davies M., Hurley J., 

Gilmore G., 1998, MNRAS, 857, 862 
Fioc M., Rocca-Volmerange B., 1997, A&A, 326, 950 
Fukushige T., Heggie, D.C., 2000, MNRAS, 318, 753 
Gieles M., 2009, MNRAS, 394, 2113 
Gieles M., Lamers H. J. G. L. M., Baumgardt H., 2008, 

IAU Symp. 246, 171 
Gieles M., Athanassoula L., Portegies Zwart S.F., 2007, 

MNRAS, 376, 809 
Gieles M., Portegies Zwart S.F., Baumgardt, H. et al., 2006, 

MNRAS, 371, 793 
Gieles M., Baumgardt H., 2008, MNRAS, 389, 28 
Gieles M., Bastian N., Lamers H.J.G.L.M., Mout J., 2005, 

A&A, 441, 949 
Giersz M., Heggie D.C., 1994, MNRAS, 270, 298 
Giersz M., Heggie D.C., 1996, MNRAS, 279, 1037 
Gnedin O.Y., Ostriker J. P., 1997, ApJ, 474, 223 
Hillenbrand LA., Hartmann L.W, 1998, ApJ, 492, 540 
Hu Y., Liu Q., Deng L., de Grijs R., 2006, in A. Vazdekis, 

Peletier R.F., eds, IAU Symp. 241, p.347 
Hurley J.R., Pols O.R., Tout C.A., 2000, MNRAS, 315, 543 
Hut P., McMillan S., Goodman J., et al., 1992, PASP, 104, 

981 
King I., 1958, AJ, 63, 109 
Kroupa P., 2001, MNRAS, 322, 231 
Kruijssen J.M.D., Lamers H.J.G.L.M., 2008, A&A, 490, 

151 



© 2009 RAS, MNRAS 000, ??-[25l 



22 Henny J.G.L.M. Lamers, Holger Baumgardt and Mark Gieles 



Kruijssen J.M.D., 2009, A&A, 507, 1409 
Kruijssen J.M.D., Mieske, S., 2009, A&A, 500, 785 
Kupper H.W., Kroupa P., Baumgardt, H., 2008, MNRAS, 

389, 889 
Lamers H. J.G.L.M., Gieles M., Bastian N. et al., 2005, 

A&A, 441, 117 
Lamers H. J.G.L.M., Gieles M., 2006, A&A Letters, 455, 

L17 
Larsen S., 2004, A&A, 416, 537 
Lee H.M., 2002, IAU Symp. 207, 584 
Leitherer C, Schaerer D, Goldader J.D., et al., 1999, ApJS, 

123, 3 
Makino J., Fukushige T., Koga M., Namura K., 2003, Publ. 

Astron. Soc. Japan, 55, 1163 
Maraston C, 2005, MNRAS, 362, 799 
McGrady N., Graham, J.R., Vacca, W.R., 2005, ApJ, 621, 

278 
Ostriker J.P., Spitzer L.Jr., Chevalier R.A., 1972, ApJ, 176, 

L51 
Portegies Zwart S.F., Hut P., Makino J., McMillan S.L.W., 

1998, A&A 337, 363 
Scheepmaker R.A., Haas M.R., Gieles M., et al., 2007, 

A&A, 469, 925 
Sommariva V., Piotto G., Rejkuba M., Bedin L.R., Heggie 

D.C., Milone A., Mathieu R.D., Moretti A., 2009, A&A, 

493, 947 
Spitzer, L. 1940, MNRAS, 100, 396 
Spitzer L.Jr., 1958, ApJ, 127, 544 

Spitzer L.Jr., 1987, Dynamical evolution of globular clus- 
ters, Princeton University Press, Princeton, N.J. 
Trenti, M., Vesperini, E. Pasquato, M., 2010, ApJ, 708, 

1598 
Vesperini E., Heggie, D.C., 1997, MNRAS 298, 898 
Von Hoerner S., 1960, ZsfAp, 50, 184 

Wielen R., 1985, in Dynamics of star clusters, IAU Sym- 
posium 113, p. 449 



are the Galactic radius, the rotation velocity and the mass 
within i?Gai- The last equality is only valid for galaxies with 
a constant rotation velocity. For clusters in elliptical orbits 
RgbI is the perigalactic distance. The half mass radius r^ of 
Roche- lobe filling clusters is 0.187 r t if the initial concentra- 
tion factor of the King density profile is Wo — 5 and 0.116 n 
if W = 7. 

It is important to estimate some of the relevant initial 
time scales of the cluster. The half mass relaxation time, t r h 



trh — 



0.138 



-rh 



3/2 



(-) 



1/2 



(lnA)" 



(A2) 



with JV ; = Mi/mi and A ~ O.llAi (Giersz & Heggie 1994), 
where mi is the initial mean stellar mass. The gravitational 
constant is G n = 0.0044985 pc 3 M© _1 Myr -2 . The crossing 
times at r^ and r t are 

^8r h 3 / 2 



t CI {r h ) = 



V0.5GM 



and 



£cr(Vt) = 



S- 



rt 



3/2 



(A3) 



(A4) 



VGMl 

The core collapse time is typically of the order of about 20 
time-dependent relaxation times. It can also be expressed in 
terms of the initial relaxation time. The BM03 models show 
that for Roche-lobe filling clusters 

0.872 



16.90 i rh L 



(A5) 



For clusters in elliptical orbits with e ^ 0.2 this expression 
underestimates the value of t cc by 20 % if the value of £ r h at 
perigalacticon is used. 

One also needs an initial rough estimate of the total 
lifetime of the cluster because dissolution depends on the 
number of stars in the cluster N — M/rn and rn depends on 
the age. The lifetime i x % can be estimated from Mi and to 
using Eq. [31] with an accuracy of about 10 percent 



APPENDIX A: RECIPE FOR PREDICTING 
THE MASS EVOLUTION OF CLUSTERS IN 
DIFFERENT ENVIRONMENTS 

Al Roche-lobe filling clusters 



a. Initial conditions 

The cluster is defined by its initial mass, Mi, initial half 
mass radius rh, initial concentration factor Wo, metallicity, 
Z, and stellar IMF. The environment is defined by only one 
parameter, to for a single cluster or £^ f for a cluster ensem- 
ble, which describes the strength of the dissolution processes 
due to the environment (see below) . Since the evolution de- 
pends on the initial cluster parameters, we first describe 
these. 

The tidal radius, r t , is 



/ M; \ 1/3 _ / GMj \ 1/3 
V2M Ga J ~ V2v Gal 2 J 



n = Rgsi x ! — — j = j — — ) i? G ai /J (Al) 
where Mi is the initial cluster mass and i?Gai , f Gal and MGai 



ti% - 3.30 x {toM - 65 } 0864 if Wo = 5 
6.27 x {toM - 80 } ' 778 if W = 7 



(A6) 



The values of to depend on the mean stellar mass, rn, which 
in turn depends on toMj 7 . So in fact t\% depends on Mi 0,61 
if Wo = 5 and Af; 067 if W = 7 (Sect.©. 

b. The total mass loss 

To calculate the mass history one has to start with mass Mi 
and then integrate numerically the mass loss rate 



^ - (iir) cv + /^ x x/ dcl a y Wx(^) cv 

+/ dclay (t) x (^M) d;s (A7) 

where the first term is the direct mass loss by stellar evolu- 
tion, the second term is the loss of stars induced by stellar 
evolution, and the last term is the dynamical mass loss (dis- 
solution) due to the tidal field, corrected for the delay in 
getting started. 
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c. The evolutionary mass loss 

The evolutionary mass loss can be calculated by means of 



m 

with 



M lum (i) *fr 



/icv(t) = 0,0 + CllX + CL2X + CI3X 



(A8) 



(A9) 



(ii) If shocks due to encounters with GMCs or spiral arms 
are the dominant dissolution process, then the value of to 
can be derived based on the properties of the spiral arms 
and the GMCs (see Gieles et al. 2006 and 2007). 
(iii) If tidal dissolution in a galaxy with a constant rotation 
velocity is the dominant effect, then to can be derived in the 
following way. 



with x = log(i/Myr). The values of the coefficients are given 
in Appendix B (Table IB2[) for a range of metallicities and 
stellar IMF^j. The luminous mass Mi um that appears in 
Eq. IA8l is the total mass M(t) minus the mass in remnants 
Afremn- The mass in remnants can be derived from the data 
in Appendix B by integrating Af rcmn from at t — 0, using 



to — t r . 



( Rgz1 \ n , /«Gal\ 1 

UxJ (1 - e) v22c7j 



(A13) 



with R in kpc and «Gai in km s _1 and t rcf = 13.3 Myr for 
Wo = 5 models and t^ ci = 3.5 Myr for Wo = 7 models. 
The mean stellar mass in the pre-core collapse phase can be 
expressed as 



^L™ - M lum (f)x[(l-/ k B «) 



d/^ E 



dt 



m = ax (t MC) 



1\b 



(A14) 



+ 



/-, rNSN d MNS ,, fWDi^WD, 

I 1 — /kick J 77 r \ L ~ ./kick J 



(If 



dt 



(A10) 



with d/i/dt > and using the resulting value of M TeulI1 (t) 
to find Mi um (i) = M(t) — Af rcmn (t). The parameters for 
Mrcmn are listed in Table IB2I for various metallicities. The 
parameter /kick is the kick factor of the remnants. If all 
neutron stars and black holes are kicked out of the cluster 
then f™ k = f ki %. — 1. In the recipe described here the kick 
factor can be adopted as a free parameter for the calculation 
of mass history of clusters. (The BM03 models described 



with (a, b, 7) = (1.528,-0.121,0.65) for Wo = 5 mod- 
els with a Kroupa IMF between 0.10 and 15 M Q , and 
(1.230,-0.094,0.80) for Wo = 7 models. (Sect. E£2j). Com- 
bining Eqs. IA13I and IA14I yields an explicit expression for 

to 



to 



prcdl 



-79 M-~^~ ,2bq 



(A15) 



above have f^ k — f ki %. = f™^ = 0, whereas the Roche- lobe with to and t Ici in Myr, M\ in M and q — (1 + 76) • For 

clusters with i prcdl M 7 < 3 Gyr the value of to is slightly 
smaller 



underfilling models have /™ = /*« = 0.90 and /g£ = 0.) 



d. The evolution-induced mass loss 

The delay function for the evolution-induced mass loss is 
/delay = 1 — exp(— t/tdciay) where the delay time scale is 

tdclay ~ 3.0 X tcr(rt). 

The scaling factor f^ x for Roche-lobe filling clusters 
can be estimated by 

5 



n 



{-0.86 + 0.40 x \og(t M?)} x (1 



(All) 



with a minimum of /j" 



0. 



e. Dynamical mass loss: dissolution 

The dissolution before core collapse can be written as 



fdMV 
V dt J 



M(t) 



1-7 



to 



(A12) 



with 7 = 0.65 or 0.80 for clusters with Wo — 5 or 7 re- 
spectively. The start of the dissolution is described by the 
function /delay (t) in Eq. lA7l 

The value of to describes the strength of the tidal field 
and other dissolution processes. So it depends on the envi- 
ronment. 

(i) If the main dissolution process is unknown, to is a free 
parameter. 

5 By applying Eq. IA9I we ignore the change in the stellar mass 
function due to the preferential loss of low mass stars. This is 
allowed because stellar evolution dominates the cluster mass loss 
at early ages before complete mass segregation is established at 

t~ 0.15ti% 



log(t prod2 ) = log(t prcdl )+0.25 [log(io prcdl A^)-3.50](A16) 

This is because clusters with shorter lifetimes contain high 
mass stars for a larger fraction of their lifetime, resulting 
in a shorter relaxation and dissolution time. For clusters in 
elliptical orbits with e ^ 0.2 the values of to are smaller 
than predicted by Eq. IA15I We found that for all models 
of clusters in elliptical orbits calculated by BM03, including 
the ones not shown in this paper, we can approximate 

log(t prodE ) = log(t prodl )-0.195 [log(t prodl Mi 7 )-3.60](A17) 

with to plcde slightly larger or smaller than to prodl , depending 
onthe value of t prcdl Afi 7 . 

The mass loss rate after core collapse can be expressed 
by a similar expression as before core collapse, except that 
it is a broken power law 



/dM\ post ' 
\~dTJ 



M(t) 



i-7c 



to" 

M ( t )l-7cc2 

*o cc2 



if M > 10 J 



if M < 10 J 



(A18) 



lO 090 t§ c for 



with 7 C c = 0.70, 7cc2 = 0.40 and with t cc 

continuity at M{t) = 10 3 . The value of tjf depends on the 

environment. 

(a) For clusters with tidal dissolution as the dominant 
effect the value of t™ depends on the mean stellar mass at 
core collapse, for the same reason as to before core collapse, 
but in this case 
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220 
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(A19) 



ref ,uf 



7.395 x (5 7 ) ' 65 if log(5 7 ) < -0.50 



(A25) 



with 7 cc = 0.70 and t^ f = 7.2 Myr if Wo = 5 and i^f = 6.2 
Myr if Wo = 7 and 

m cc = a (t Mr l ) b (A20) 

with (a, b) = (1.585,-0.0984) for W^o = 5 models, and 
(0.893, -0.0691) for Wo = 7 models. Clusters in elliptical 
orbits have a smaller mean stellar mass 



rn.cc = m t 



: ( e = 0)x(l-e) u 



(A21) 



because of their shorter lifetime. Substitution of Eq. lA20l or 
[MT\ into lA"T9l yields the value of io c . 

(b) For clusters in an environment where external effects 
(other than the tidal field) dominate, the value of £q c can De 
derived from to by applying the same jump in mass loss at 
tec as for dissolving clusters. This implies 

¥ ^^~ (A22) 



t() 



We now have a full set of equations that describe the 
mass loss rates due to stellar evolution, evolution-induced 
mass loss and dissolution before and after core collapse of 
Roche-lobe filling clusters. With these sets of equations the 
mass evolution of Roche-lobe filling clusters can be calcu- 
lated by means of numerical integration of Eq. IA7I 

A2 Roche-Lobe underfilling clusters 

The recipe for calculating the mass history of clusters that 
are initially Roche-lobe underfilling proceeds along the same 
lines as that for the Roche-lobe filling clusters with a few 
modifications. The cluster is defined by the same initial pa- 
rameters as in Sect. I All plus the underfilling factor 5w de- 
fined by Eq. 1231 with 5=1 for Roche-lobe filling clusters. 
Clusters with an initial density concentration of Wo =£ 7 and 
5 < 1 quickly expand and reach a density distribution close 
to Wo < 7. For that reason $7 is also needed to describe the 
mass history. 

The values of the initial parameters and time scales are 
the same as described in Sect. IA1I The calculation of the 
mass loss rate also proceeds along the same lines. The scaling 
factor /iJ'cf for the induced mass loss is 



fST = -0-86 + 0.40 x log(toM^) + 2.75 log 5s 



(A23) 



The delay time of clusters with —1.0 < log 5s < —0.20 can 
be estimated by 



-f /I Q1 m-3 v, IK \- 1-989 . 1.605 

tdelay = 4.31 10 X (#5) trf, 



(A24) 



For log ^5 > —0.20 the delay time is about the same as for 
Roche-lobe filling clusters, idday — 3.0t cr (r t ). 

The dissolution is described by Eq, IA12l but in this case 
7 = 0.80 even if the cluster started with Wo — 5. For clusters 
in an environment where tidal stripping is not the dominant 
mass loss mechanism, to is a free parameter. For clusters 
with tidal stripping in a galaxy with a constant rotation 
velocity the dissolution parameter is given by Eq. lA13l with 
& U f = 3.50 Myr if log(3V) > -0.50 and 



The mean stellar mass in the pre-core collapse phase is 
given by 

m prc -cc = a x {t M?') b x 5? (A26) 

with (a, b, c, 7)= (1.38, -0.0984,+0.101, 0.80) for a Kroupa 
mass function of ro m i n = 0.10 and m max = 100 Mq. For 
other values of m ma x the parameters will be about the same 
because the most massive stars are lost before dynamical 
effects are important. Combining Eqs, IA131 and 1X261 gives 
an explicit expression for to 



to 



"Tef.uf A 



(W" 



220 

WGal 



-11 J1J 1 bq rg--ycq 



with q = 1/(1 + 76). 

The core collapse time is given by 

t C c — OZ.U X t r h X t^5 



(A27) 



(A28) 



The mass loss rate after core collapse is given by Eq. 
IA18I The value of t(j c can be derived from t(j c = t^ f cc x 
m cc "^ cc with m cc = a cc (toMj 7 ) 6 " 87°° with (a cc , 6 CC , c cc ) = 
(1.507,-0.0984,0.207) (see Eq-HSJ and 

C™ = 7.50 x f 5 127 (A29) 

The value of 7cc = 0.70 if M > 10 3 and 7 cc2 = 0.40 if M(t) < 
10 Mq. This change is due to the variation of the Coulomb 
logarithm for small numbers of stars (see Sect. I6.2f) . The 
dissolution parameter changes from t™ to t cc2 = 10°' 90 to c 
if M(i) < 10 3 M Q (see Eq. [A"T8|) . 

With these equations the mass history of initially 
Roche-lobe underfilling clusters can be calculated. 



APPENDIX B: MASS LOSS AND REMNANT 
PRODUCTION BY STELLAR EVOLUTION 

In this section we describe power law approximations for 
calculating the different contributions to the mass of a clus- 
ter due to stellar evolution as a function of age for non- 
dissolving clusters. This is done for metallicities of Z = 
0.0004, 0.001, 0.004, 0.008 and 0.02 for a Kroupa-IMF from 
rrimin < m < m max with m max = 100 Mq and m m i n = 
0.1 Mq. The mean initial mass is < m >i= 0.638 Mq. 
We adopted the stellar evolution described by Hurley et al. 
(2000). 

Table IB1I gives the relation between age and the mass 
of the star that ends its life at that age. Table lB2l gives the 
parameters for calculating the variation of several cluster 
quantities by means of polynomial fits of the type y = 00 + 
aiX + a'zx 2 + a-iX 3 with x — log(t/Myr). These quantities 
are: p(i) = M(t)/Mi is the remaining mass fraction of the 
cluster; /ZBH(t) = Mbh/Mi, /UNs(t) and /AvD(t) are the mass 
fraction (relative to the initial cluster mass Mi) of the black 
holes, neutron stars and white dwarfs; < m >, < m >bh, 
< m >ns and < m >wd are respectively the mean masses 
of all stars and of the black holes, the neutron stars and 
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Table B 1 . The maximum initial mass of a star that has survived 

at age t is given by the polynomial: 

log(m max ) = ao + aiy + a2y 2 + asy 3 + asy 5 with y = log(t/Gyr). 



z 


tt 


0.1 


ai 


0.3 


a r> 


0.0004 


0.2732 


-0.3864 


0.05628 


0.01524 


-0.005902 


0.0010 


0.2940 


-0.3892 


0.05075 


0.02055 


-0.006962 


0.0040 


0.3078 


-0.3910 


0.04709 


0.02407 


-0.007662 


0.0080 


0.3273 


-0.3867 


0.04203 


0.02740 


-0.008361 


0.0200 


0.3463 


-0.3789 


0.03389 


0.03290 


-0.009552 


0.0500 


0.3307 


-0.3653 


0.02304 


0.03663 


-0.010346 



the white dwarfs (in units of Mq) and log(L(t) / Lq) is the 
luminosity of the cluster, normalized to an initial cluster 
mass of 1 Mq. The powerlaw fit for the luminosity is split 
in two parts for two age ranges. The resulting power law fits 
can be used to calculate the changes in the mass functions 
of dissolving clusters due to stellar evolution. 

From the data of this table one can also derive the num- 
ber of black holes, neutron stars and white dwarfs as a func- 
tion of time 

iVco(t) = Mi Hco(t)/ < m > ro (t) (Bl) 

where "co" stands for compact object, i.e. BH, NS, or WD. 
Similarly, the total number of objects, stars plus compact 
objects, is N — Mi fi/ < m >. The number of luminous 
(non-compact) stars is JVi um = N — Nbh — -/Vns — -Wwd and 
the mean mass of the luminous stars is 

< m >lum= Mi (n(t) — flBH — jUNS — MWD)/Wlum (B2) 

If compact objects are lost from the cluster by their kick 
velocities, then the values can be adjusted by applying kick- 
factors / k ™ k , f^ k and fg£ with / k ™ k = 0.90 if 90% of the 
black holes are lost from the cluster during their formation. 
For instance JVbh = (1 — fkuL) Mi /^bh/ < rn >bh- 

For clusters with a Kroupa IMF with a higher lower 
mass limit, m m i n , the parameters can be calculated by sub- 
tracting the mass of the missing lower part of the IMF from 
the normalization to Mi. This is valid because the low mass 
stars with m < O.8OM0 do not contribute to the evolu- 
tion at t < 50 Gyr. For instance, a cluster of a given initial 
mass Mi and an IMF in the range of 0.50 < m < 100 M @ 
has 1.091 times as many stars with m > O.5A^0 as a clus- 
ter with the same Mi but with an IMF in the range of 
0.5 < m < 100 Mq. So the fractions fi, fiBH, Mns and (iwd 
will be higher by 1.091. 

Similarly, for cluster with a Kroupa IMF, but with a 
lower value of m max the mass fractions fi, /^bh, Mns and 
/iwD can be derived as follows. Calculate first the age of 
the cluster, £(m max ), at the time when a star of mass m ma x 
ends its life from Table |BT1 and the value of ^(£(m max )) from 
Table lB2l At any time t subtract the value of ^BH(t(w m ax) 
from fiBH(t) and normalize it to the new initial mass by 
dividing the result by /i(£(m max )) to find the new value of 
jUngf. The same method can be used to find the values of 
/^ ncw , MneT an d Mwd- The number of compact stars is found 
by iV co = N co (t) — AT co (t(m m ax)) and the mean mass is < 
m> co =Mi/x? ew /Af C o. 
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Table B2. Parameters for calculating various stellar evolution parameters for a cluster with a Kroupa IMF in the range of 0.10 < m < 
100 M , as function of metallicity Z and age t by means of a polynomial: 
y = ao + a\x + a,2X 2 + azx s with x = log(i/Myr). 



y 


Z 


«o 


0.1 


a 2 


0.3 


a 


■*-miii 


^max 


2/l^max/ 


!-'■ 


0.0004 


1.0541 


-0.10912 


-0.01082 


0.00285 


0.000 


1.0000 


4.6990 


0.5981 


I- 1 


0.0010 


1.0469 


-0.10122 


-0.01349 


0.00306 


0.000 


1.0000 


4.6990 


0.5909 


H 


0.0040 


1.0247 


-0.08307 


-0.01845 


0.00336 


0.000 


1.0000 


4.6990 


0.5756 


I-' 


0.0080 


1.0078 


-0.07456 


-0.02002 


0.00340 


0.000 


1.0000 


4.6990 


0.5682 


A* 


0.0200 


0.9770 


-0.05709 


-0.02338 


0.00348 


0.000 


1.0000 


4.6990 


0.5536 


MBH 


0.0004 


-1.0711 


2.66839 


-2.07015 


0.52230 


0.000 


1.0000 


1.1406 


0.0543 


Mbh 


0.0010 


-1.1254 


2.83983 


-2.25098 


0.58431 


0.000 


1.0000 


1.1141 


0.0525 


Mbh 


0.0040 


-0.7770 


2.00007 


-1.60085 


0.42042 


0.000 


1.0000 


1.1388 


0.0455 


/^BH 


0.0080 


-1.2000 


3.26274 


-2.86105 


0.83466 


0.000 


1.0000 


1.1208 


0.0380 


Mbh 


0.0200 


0.0274 


0.00000 


0.00000 


0.00000 


0.000 


1.0000 


1.1000 


0.0274 


Mns 


0.0004 


-0.1343 


0.25061 


-0.15593 


0.03497 


0.020 


1.1226 


1.7885 


0.0152 


A«NS 


0.0010 


-0.1099 


0.19575 


-0.11522 


0.02497 


0.020 


1.1236 


1.7803 


0.0143 


MNS 


0.0040 


-0.0849 


0.14542 


-0.08131 


0.01723 


0.000 


1.1136 


1.7120 


0.0122 


MNS 


0.0080 


-0.0367 


0.04471 


-0.01078 


0.00071 


0.022 


1.0840 


1.7113 


0.0118 


A'NS 


0.0200 


-0.0128 


-0.00396 


0.02480 


-0.00823 


0.001 


1.0091 


1.6395 


0.0111 


CWD 


0.0004 


0.0064 


-0.03566 


0.01882 


-0.00034 


0.122 


1.7571 


4.6990 


0.2191 


CWD 


0.0010 


0.0087 


-0.03363 


0.01680 


-0.00007 


0.110 


1.7120 


4.6990 


0.2144 


CWD 


0.0040 


0.0014 


-0.01864 


0.01000 


0.00056 


0.056 


1.6300 


4.6990 


0.1927 


CWD 


0.0080 


0.0031 


-0.01809 


0.00916 


0.00058 


0.048 


1.6020 


4.6990 


0.1805 


CWD 


0.0200 


-0.0173 


0.00709 


0.00025 


0.00133 


0.012 


1.5928 


4.6990 


0.1595 


< m > 


0.0004 


0.6681 


-0.06910 


-0.00686 


0.00180 


0.000 


1.0000 


4.6990 


0.3787 


< m > 


0.0010 


0.6635 


-0.06416 


-0.00855 


0.00194 


0.000 


1.0000 


4.6990 


0.3745 


< rn > 


0.0040 


0.6495 


-0.05264 


-0.01170 


0.00213 


0.000 


1.0000 


4.6990 


0.3648 


< m > 


0.0080 


0.6388 


-0.04727 


-0.01269 


0.00216 


0.000 


1.0000 


4.6990 


0.3606 


< rra > 


0.0200 


0.6193 


-0.03621 


-0.01481 


0.00220 


0.000 


1.0000 


4.6990 


0.3504 


< rn >bh 


0.0004 


33.5859 


-16.14886 


0.00000 


0.00000 


0.000 


1.0000 


1.0925 


15.9433 


< ™ >bh 


0.0010 


31.5692 


-14.57680 


0.00000 


0.00000 


0.000 


1.0000 


1.0928 


15.6397 


< rn >bh 


0.0040 


30.2147 


-15.18139 


0.00000 


0.00000 


0.000 


1.0000 


1.0770 


13.8643 


< m >bh 


0.0080 


23.6200 


-11.10258 


0.00000 


0.00000 


0.000 


1.0000 


1.0647 


11.7992 


< rn >bh 


0.0200 


9.3469 


0.00000 


0.00000 


0.00000 


0.000 


1.0000 


1.0500 


9.3469 


< m > NS 


0.0004 


16.3368 


-26.6599 


16.46480 


-3.46414 


0.000 


1.1121 


1.7650 


1.5265 


< "i > NS 


0.0010 


13.4717 


-20.6296 


12.27190 


-2.50150 


0.000 


1.1121 


1.7640 


1.5367 


< m > NS 


0.0040 


9.8973 


-13.4125 


7.34369 


-1.37739 


0.000 


1.0747 


1.7315 


1.5403 


< m > NS 


0.0080 


14.1365 


-22.1270 


13.19290 


-2.66813 


0.000 


1.0747 


1.7211 


1.5308 


< m > NS 


0.0200 


13.4243 


-20.6297 


12.06450 


-2.37325 


0.000 


1.0500 


1.7020 


1.5601 


< rn >wd 


0.0004 


2.8955 


-1.21969 


0.25054 


-0.01995 


0.001 


1.7847 


4.6990 


0.6263 


< m >wd 


0.0010 


3.0223 


-1.33991 


0.28235 


-0.02259 


0.001 


1.7847 


4.6990 


0.6167 


< m > WD 


0.0040 


3.0228 


-1.34601 


0.27457 


-0.02079 


0.000 


1.7473 


4.6990 


0.6035 


< m >wd 


0.0080 


3.0213 


-1.38321 


0.29001 


-0.02239 


0.000 


1.7099 


4.6990 


0.6021 


< m > WD 


0.0200 


2.8737 


-1.25599 


0.24619 


-0.01755 


0.000 


1.6352 


4.6990 


0.5869 


log L/Lq 


0.0004 


4.3507 


-2.72949 


0.88598 


-0.14060 


0.000 


1.0000 


2.5272 


0.8419 


log L/L 


0.0010 


4.4074 


-2.78762 


0.89986 


-0.14251 


0.000 


1.0000 


2.2258 


1.0893 


log L/L 


0.0040 


4.9761 


-3.84505 


1.49603 


-0.25205 


0.000 


1.0000 


1.8173 


1.4165 


log L/Lq 


0.0080 


4.5921 


-3.19055 


1.11463 


-0.18185 


0.000 


1.0000 


2.7372 


0.4807 


log L/L 


0.0200 


4.8125 


-3.62327 


1.30573 


-0.20546 


0.000 


1.0000 


2.8070 


0.3860 


log L/L GJ 


0.0004 


3.2921 


-1.20074 


0.13311 


-0.01647 


0.000 


2.5272 


4.6990 


-1.1199 


log L/L 


0.0010 


5.4132 


-3.20696 


0.73195 


-0.07364 


0.000 


2.2258 


4.6990 


-1.1351 


log L/L 


0.0040 


4.1575 


-2.04685 


0.35860 


-0.03425 


0.000 


1.8173 


4.6990 


-1.0962 


log L/Lq 


0.0080 


2.8690 


-1.35525 


0.27530 


-0.03615 


0.000 


2.7372 


4.6990 


-1.1713 


log L/L 


0.0200 


10.4008 


-7.65806 


1.98530 


-0.18815 


0.000 


2.8070 


4.6990 


-1.2696 



(a) cr is the standard devaiation of the fit in the range of x m [ n < x < x max . 

(b) x m i n and x max indicate the validity range of the approximation. 

(c) All values of ii are expressed in units of M\, so ii = M(t)/M[, /ibh = ^BH/A^i etc. 

(d) The values of /ibh and /^ns an d the values of < m >bh an d < rn >ns are constant fot axi.x mv c^t.u 

(e) < m > is the mean stellar mass, < m >bh is the mean stellar mass of the black holes; etc, in umfe o 

(i) The values of loslL/Lr.-)) are for a cluster of M; = IMo. The fit is split in two age ranges. 
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